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Synopsis 


A number of real-life systems, like a band-saw, a magnetic tape, a power transmission 
belt, an overhead cable for aerial tramways, a flexible manipulator arm having prismatic 
joint etc., can be modelled as axially travelling continua. With increasing axial speed, 
the productivity increases but the performance deteriorates due to the transverse vibra- 
tion caused by external or parametric excitations. The separation-of-variable method, 
a well-known tool for the analysis of an axially stationary member, cannot be applied 
for a travelling member due to the presence of the ‘gyroscopic’ term. Further, unlike in 
a stationary medium, the natures , of the propagating transverse waves in a travelling 
medium become complicated to form a clear physical idea. These complications are 
further enhanced in a dispersive medium. 

In addition to the above difficulties, the analysis also becomes cumbersome due 
to the presence of geometrical non-linear terms, which are significant during large 
oscillations at high axial speeds. Such large oscillations are caused by resonance or hard 
non-resonant excitations. The non-linear terms are important not only for obtaining 
quantitatively accurate resuls but also for explaining several phenomena which cannot 
be predicted from a linear analysis. Furthermore, any measure based upon a linear 
model to control the unwanted vibration may turn out to be ineffective in the presence 
of non-linear terms. 

The aim of the thesis is to study the effects of non-linearities on the transverse 
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vibration of an Euler-Bernoulli beam moving over two simple supports. Various op- 
erating conditions generating different kinds of excitations are considered. Compared 
to the vast amount of literature available on the linear vibration of travelling continua 
(discussed in Chapter 1), the non-linear oscilations of a travelling beam have so far 
received less attention. 

The non-linear equation of motion is derived in Chapter 2. Since the small non- 
linear term in the equation of motion is considered as a perturbation to the linear 
equation, a detailed knowledge of the linear system is essential. The rest of Chapter 
2 is devoted to study the free and forced linear vibrations by two apparently distinct 
methods, namely the ‘modal’ analysis and the ‘wave-propagation method’. While 
the discussion of the former is limited to the existing results, the latter method is 
developed for the first time. The numerical results obtained from these two methods 
are compared. Damping is not considered unless the presence of the damping term 
changes the qualitative behaviour of the system. 

In Chapter 3, a simple and computationally efficient method of obtaining the near- 
resonance response of a non-linear continuous system is presented. The method is 
elucidated, for the sake of simplicity, with an axially stationary beam having either 
simply-supported or clamped-clamped boundary conditions. The non-linear normal 
modes of the beam are derived first by perturbing the linear normal modes and then also 
using the wave-propagation method. These non-linear normal modes are subsequently 
used to obtain the near-resonance response. The results so obtained are then compared 
with the experimental results available in the literature. The efficacy of the present 
method over the Galerkin’s technique is highlighted. 

The results of Chapters 2 and 3 are extended in Chapter 4 to derive the response of 
the non-linear travelling beam under both resonant and hard non-resonant excitations. 
A concept of non-linear complex normal mode is first presented. The near-resonance 
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response is subsequently obtained using these non-linear complex normal modes. The 
methods of deriving both the non-linear complex modes and the near-resonance re- 
sponse closely resemble those outlined in Chapter 3. Like in any other non-linear sys- 
tem, multiple values are obtained for the near-resonance response. However, all of these 
values are not physically realizable. The physically observable results are delineated 
by carrying out a stability analysis. When the excitation is hard non-resonant one, the 
modal method encounters difficulties. The method of wave-propagation, derived for 
the linear system in Chapter 2, has been effectively extended for such a situation. 

Chapter 5 deals with the effects of the non-linear term on a parametrically ex- 
cited travelling beam. The source of the parametric excitation is taken to be the 
tension fluctuation due to the unbalance in one of the pulleys which is mounted on an 
axially-flexible support. Single or multiple limit cycles, generated by the combination 
of non-linearity and parametric excitation in the presence of small viscous damping, are 
obtained using the non-linear complex normal modes. The stable limit cycles are iden- 
tified by carrying out a stability analysis. The response of the beam to simultaneous 
parametric and external harmonic excitations is also obtained. The non-linear complex 
normal modes are used for this purpose. Again a stability analysis for determining the 
physically observable solutions is presented. 

The stability of an accelerating or decelerating non-linear beam is presented in 
Chapter 6. This situation arises during the starting or stopping phases. A Lyapunov’s 
stability analysis is first carried out to study the stability characteristics of the beam 
having any arbitrary acceleration. For a decelerating beam, the stability is determined 
using a perturbation technique based on Multiple Time Scale (MTS) method. The 
uniform deceleration which alters the magnitude of the uniform axial speed but not 
the direction of axial movement is considered as a small perturbation. 

A controller to reduce the vibration of a non-linear beam moving with uniform 
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axial speed is discussed in Chapter 7. A passive controller in the form of a roller guide 
is proposed. The finite compliance of the guide in the transverse direction, and the 
frictional force appearing at the guide-beam interface and acting in the axial direction 
are incorporated. The results of the linear analysis, for both free and forced vibra- 
tions, are used to obtain the non-linear complex normal modes that are subsequently 
utilised to get the near-resonance response. The effectiveness of the roller in reducing 
the amplitudes of both linear and non-linear vibrations of the beam under harmonic 
excitation has been brought out with the help of numerical examples. 

The main conclusions of the present work are summarized in Chapter 8. Some 
directions for further research are also mentioned therein. The major conclusions of 
the analyses presented in the thesis can be summarized as listed below: 

(i) The natures of the waves associated with the free modal vibration of a travelling 
beam change depending upon the relations between various system parameters. During 
linear modal vibration, either two propagating and two evanescent waves (one of each 
kind in both upstream and downstream directions) or four propagating waves (one in 
the upstream direction and the remaining three in the downstream direction) appear 
depending on the axial speed. The well-known ‘phase-closure’ principle is satisfied for 
both the cases. The ‘phase-closure’ principle is also satisfied during non-linear modal 
oscillations of such a beam. 

(ii) The near-resonance response to an external harmonic and/or parametric excitation 
can be obtained easily using the non-linear complex normal modes. This method is 
computationally more efficient than the commonly-used Galerkin’s technique. 

(iii) The steady-state response of the beam to a non-resonant hard excitation can be 
obtained by wave-propagation analysis. In this method, a closed form transfer function, 
i.e., the Laplace transform of a point-impulse response, is first obtained. The linear 
response obtained using this transfer function is then perturbed to get the effects of 
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the small non-linear term. The present method is shown to be more accurate and 
computationally more efficient than either Galerkin’s technique or non-linear complex 
modal analysis. 

(iv) The non-linear natural frequencies of a travelling beam increase with the amplitude 
of vibration. The ‘jump’ of the amplitude of the response envelope can be observed un- 
der a near-resonant excitation. Further, one or multiple ‘limit cycle’s of the amplitude 
of the response envelope are possible when the beam is parametrically excited. These 
three phenomena show striking resemblance with those exhibited by a hard Duffing 
oscillator. 

(v) For a parametrically excited travelling beam with an external harmonic excitation, 
the effects of various system parameters depend on the relative strengths of the two 
forms of excitation. Moreover, a particular phase relationship between the two forms 
of excitations (whose frequencies are related) results in a minimum value of the steady- 
state response. 

(vi) A continuously accelerating beam remains stable, but for a decelerating beam in- 
stability may appear depending upon the magnitudes of damping and deceleration. 
However, the amplitude of a decelerating beam does not grow unboundedly and the 
beam regains its stability in the long run. The non-linear term does not play any role 
so far as the stability is concerned. It merely changes the frequencies of oscillations. 

(vii) An intermediate guide can be effectively used as a passive controller of vibration. 

4 

The choice of the suitable guide-location plays a very important role in reducing the 
level of vibration. The guide also has a stabilizing effect so far as the divergence in- 
stability is concerned. But the frictional force between the guide and beam may add 
to instability. Thus, to use the controller effectively, the friction has to be minimised. 
Towards this end, a roller-guide is suggested. 
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Chapter 1 


INTRODUCTION 


1.1 Introduction 

A number of appliances in diverse fields of engineering are modelled as axially mov- 
ing continua. Some of the common examples are saw-bands used in forest industries, 
magnetic tapes in tape-recorders and computers, travelling threadlines in textile in- 
dustries, belt-drives in automobiles, flexible manipulator arms having prismatic joints, 
etc. High productivity often requires these systems to operate at the highest possible 
speed. However, at a high speed, the performance deteriorates primarily due to vibra- 
tions of the moving element in a direction orthogonal to that of the axial movement. 
The presence of the ‘gyroscopic term’ in the equation of motion of transverse oscillation 
of such systems makes the problem quite different from that of an axially stationary 
slender member. For the latter case, a great variety of mathematical tools have been 
developed over the centuries. Further, the close relationship between the vibration 
of and wave propagations in axially stationary members has helped the researchers 
to gain insight into their oscillatory behaviour. But even a simplistic linear model of 
a travelling continuous system presents some difficulties as the separation of variable 
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method, a well-known tool for solving a boundary value problem, is not applicable in 
this situation. The natur&-©f the propagating transverse waves in a travelling medium, 
especially if it is dispersive in nature like a beam, becomes complicated to form a clear 
physical idea. These difficulties have motivated researchers to give special attention to 
this class of systems. 

Apart from the technical difficulties mentioned above, the real situation is further 
complicated by the inevitable presence of various non-linear terms at different levels of 
approximation. Depending upon their sources, these non-linear terms are of two kinds:, 
one appearing due to material non-linearities, i.e., non-linear variation of stress with 
strain, and the other resulting from geometric non-linearities i.e., non-linear strain- 
displacement relationship. While the former depends upon the material properties of 
the vibrating medium, the latter is encountered during large-amplitude oscillations, for 
example, due to resonance or non-resonant hard excitations. The importance of the 
geometrical non-linearities increases with increasing axial speed of a travelling member. 
The results of a non-linear analysis not only quantitatively differ : from that obtained by 
a linear analysis, certain phenomena observable in real life can never be explained from 
a linear analysis. Furthermore, a control measure based only upon the approximate 
linear model may turn out to be ineffective in the presence of non-linear terms. 

For a continuous system, the non-linear analysis becomes increasingly difficult due 
to modal coupling. The linear vibration analysis of a continuous system is usually 
carried out by mathematically discretizing the equation of motion, to a set of infinitely 
many decoupled simple equations, with the help of orthogonal normal modes. Such 
modal decoupling is generally not possible in the non-linear case. Some simplifications 
are achieved during resonant vibrations, when the contributions to the overall response 
from one or few linear modes are orders of magnitude higher than those from the 
others. But the analysis becomes cumbersome during a non-resonant large oscillation. 
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especially due to a high-frequency excitation, when a good number of modes contribute 
significantly to the response. 

The aim of this thesis is to study the effects of non-linear terms on the vibrations 
of a travelling beam operating under various conditions. To this end, special analytical 
tools for solving the problems under consideration have been devised. The physical 
understanding of the oscillatory behaviour has been enhanced by explaining both the 
free and forced vibrations in terms of various waves propagating in this dispersive 
medium. An easily implementable passive controller has also been suggested to reduce 
the non-linear vibration of such a system. 


1.2 Literature Review 

Following Aitken’s earliest experiment on the vibrations of a travelling belt, researchers 
have been trying to understand and control the dynamics of travelling continuous 
systems [1]. These efforts have been compiled in several review articles [2]- [6]. A 
very rich literature also exists for the closely related problem of vibration of pipes, 
where the axial material movement is maintained by the fluid flow [7]. 

As mentioned in Section 1.1, the vibration of a moving medium, whether supported 
at one or both ends, differs from that of an axially stationary medium. While both the 
mass and energy of the span are conserved in a stationary member, the energy of a 
travelling member is not conserved [8]- [11]. Further, if the support is only at one end, 
the length of the span and hence, the oscillatory mass changes as the member moves. 
Thus, the axially moving systems can be divided into two categories, one having ‘fixed 
span’ and the other with ‘variable span’. These two kinds of system have been treated 
separately in the literature. 
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1.2.1 Free Vibration 


1. 2.1.1 Strings and Cables 

Free linear vibration of a string travelling over two supports has been studied analyti- 
cally [12]- [16]. The theory of wave propagation has also been used [9, 17]. Following 
the works of Meirovitch [18, 19] and other researchers [20] in the field of linear gyro- 
scopic systems, the modal oscillation of a travelling string has been studied in terms of 
the complex normal modes. Here the waves in the upstream and downstream directions 
travel with different speeds. Consequently a standing envelope rather than a standing 
mode is formed. The impossibility of the standing modes is also evident from the 
equation of motion, where the temporal and the spatial variables are coupled through 
the Coriolis component in the transverse force, known as the ‘gyroscopic’ term in the 
literature. The centrifugal force, on the other hand, compels the natural frequencies to 
decrease with increasing axial speed until a ‘critical’ speed is reached when all of them 
vanish simultaneously. At this critical speed, no wave can propagate in the upstream 
direction. 

Systems like continuously-supported conveyer belts, air-guided magnetic tapes, 
translating paper-pulp sheets supported by air-jets etc. are modelled as an axially 
moving string having intermediate spring supports. The nature of the propagating 
waves and thus the d 3 mamics of such systems are quite different from those of the 
systems supported only at the ends of the span [21]- [27]. For example, in a travel- 
ling string elastically supported over a part of the span, the harmonic nature of the 
time-response is lost if the frequencies of the propagating waves become smaller than 
a ‘cut-off’ frequency [21, 28]. However, the ‘critical’ speeds of these systems are not 
altered. On the other hand, the dry frictional force present in the intermediate sup- 
ports like capstans, cylinders, eyelets, read/erase heads, rollers etc., reduces the critical 
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speed significantly [29]. Studies with various other constraints reveal that while a pure 
elastic constraint increases the natural frequencies, an inertia constraint in the form 
of a stationary point mass has a reverse effect [26, 30]. The well-known ‘eigenvalue- 
inclusion principle’ holds good for both these cases [30]. The complicated situation 
with both kinds of constraints has also been analysed [31]. 

The dynamics of a translating cable has also received attention [32]- [36]. Due to the 
presence of a small sag, the natural frequencies of such a system do not always decrease 
monotonically with increasing axial speed. The angle of inclination of the cable, caused 
by a difference in the heights of the end supports, may change the dynamics significantly 
[33]. For a small axial speed, the equilibrium configuration of a sagged cable is unique 
and is known as the ‘minimum catenary’. But for a sufficiently large translating speed, a 
second arch-like equilibrium configuration, known as the ‘maximum catenary’, becomes 
stable [35, 36]. 

1.2. 1.2 Beams 

The vibration problem of a travelling beam was first studied in 1965 with reference 
to a saw-band which was modelled as an Euler-Bernoulli beam travelling over two 
frictionless simple supports [37]. This and subsequent analyses have shown that the 
natural frequencies of such a beam, like those of a string, decrease with increasing axial 
speed [38]- [40]. But, unlike in a string, there exist different critical speeds for different 
complex modes of the travelling beam. The compressive load, generated due to the 
centrifugal force at any of these critical speeds, becomes equal to the Euler buckling 
load for the corresponding mode. These critical speeds, however, increase if the beam is 
driven by pulleys having finite support compliance in the axial direction [37]. Studies 
on the vibrations of travelling beams having non-trivial equilibrium configuration have 
also been carried out [41, 42]. The static deflection of the beam from the trivial con- 
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figuration is due to manufacturing flaws, finite end-curvatures because of the pulleys, 
etc. The natural frequencies of the odd order modes of such a beam do not vanish at 
any speed, though the natural frequencies of the even order modes disappear at the 
critical speeds. 

For an unguided saw-band driven by pulleys, the vibrational energy does no longer 
remain confined to any of the two identical spans but gets transferred from one span 
to the other periodically [43, 44]. Consequently, introduction of damping in one span 
helps to damp out the oscillations in the other [45]. The periodic transfer of energy 
among the spans, known as the ‘beating phenomenon’, is absent if the band is replaced 
by a srting [43]. The band/wheel system has been studied using two different models. 
In one model, the non-trivial static equilibrium configurations of the spans, due to the 
finite end-curvatures at and near the pulleys, are taken into account and the linearized 
equations of motion of both the transverse and longitudinal vibrations of the two spans 
have been obtained [44]. The longitudinal vibration, which affects the pulleys, plays 
a significant role in coupling the oscillatory motions of the two spans. The other 
simplified model assumes the trivial equilibrium configurations of the spans and. the 
coupling is manifested through the slope- and monient-continuity at each adj ascent 
end [46, 47]. It may be mentioned that in reference [46], the boundary conditions 
are incomplete ( only six of the required eight conditions are presented). Analyses of 
both the above models show that the frequencies of a band/wheel system appear in 
groups, each consisting of two closely separated frequencies. The difference between the 
two frequencies of any such group governs the frequency at which energy is exchanged 
between the spans. The effects of various parameters like band tension, axial speed, 
pulley radii, etc., on the vibration coupling have been reported. The coupling decreases 
with increasing pulley radii and/or the tensions in the bands. As the axial speed is 
increased, the coupling first decreases until the spans get decoupled for a particular 
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mode (i.e., two frequencies belonging to a group become identical). If the axial speed 
is increased further the coupling gets stronger. The decoupling corresponds to the 
intersection of the frequency vs. axial speed plots. However, a slight disorder in the 
spans, for example, due to a difference in the initial band tensions, may introduce a 
new phenomenon, known as ‘frequency-loci veering’ or ‘curve veering’ [48]. As the 
name suggests, the two frequencies, in stead of becoming equal, abruptly veer away 
from each other. The aforesaid curve veering results in ‘vibration localization’ or ‘mode 
localization’ in one of the two spans [47, 49]. This kind of localization of oscillatory 
motion has been observed in various physical and engineering systems consisting of 
periodic elements with disorder [50]- [54]. It may be worthwhile to note that the 
localization occurs also in a friction guided travelling string [29]. 

1.2. 1.3 Strings/Beams with Attached Systems 

The dynamics of a travelling member changes significantly if a rigid or fiexible body 
is attached thereto. Examples of such a system include an aerial tramway, a ski-lift, 
a monocable ropeway, a conveyer belt, a saw-band having material imperfections like 
weld-melts, etc. The non-harmonic nature of the response during free vibration of such 
a system has been confirmed by both exact and approximate analyses [55]- [57]. In 
the exact analysis, the unknown force between the string and the appendage is solved 
by formulating a Volterra equation of the first kind. In the approximate method, the 
inhomogeneity in the mass distribution is considered as a small perturbation to the 
original system. Theoretical as well as experimental results show that the amplitude 
of vibration of a string having an attached body is more than that of a string without 
any appendage. Thus, for safe operation of such a system, the axial speed must be 
restricted to a value much lower than the critical speed. 
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1.2.2 Vibration due to External/Parametric Excitations 

The responses of a constrained or unconstrained travelling continuous system to 
force/displacement excitations have been obtained by three different methods, namely, 
Galerkin’s technique or the assumed mode method [58], the exact modal analysis 
[16, 59] and the transfer-function method [60]- [63]. In the first method, the equations 
of motion are discretized, with the help of the modes of an axially stationary member, 
to a set of ordinary differential equations in temporal coordinates which are then solved 
to get the response. However, the implicit assumption of the separation of the spa- 
tial and temporal variables results in slow convergence of the series solution. A faster 
convergence can be achieved by discretizing the equation of motion to several complex 
ordinary differential equations with the help of the exact complex normal modes [64]. 
The series form of the response can be altogether avoided by formulating closed-form 
transfer function (i.e., the Laplace transform of the Green functions) of the travelling 
member [60, 62]. The response to any arbitrary excitation is then derived with the 
help of this transfer function. The last method enjoys an added capability of tack- 
ling any arbitrary boundary condition. For a travelling string, different terms of the 
transfer function have been also explained in terms of various propagating waves [63]. 
The response of the string to a displacement excitation can be very easily obtained 
using the wave-propagation theory. For a beam, however, the waves are significantly 
distorted making the wave-propagation analysis difficult. It should be mentioned here 
that if the load is time-invariant, the problem of wave propagation does not arise and 
the exact response can be easily obtained [65]. 

Apart from direct external excitations, the travelling members are often subjected 
to parametric excitations. Such excitations are caused by the fluctuations in either 
tension [66, 67] or axial speed [68]- [71] of the members. For a travelling string 
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(especially V-belts), the tension fluctuations are caused mainly by the belt-defects, 
whereas the eccentricities of one or both the pulleys may periodically change the tension 
of a saw-band. The parametric fluctuations in the belt-drive of an automobile occur 
through velocity modulations due to the fluctuating turning moment of the engine. 
Various aspects of a parametrically excited travelling system have been reported in 
the literature [66]- [71]. Both primary and secondary resonances have been studied. 
The response of a travelling beam having random tension fluctuations has also been 
obtained using stochastic averaging technique [72]- [74]. 

1.2.3 Control 

It is well known that the vibration deteriorates the overall performance of travelling 
systems. Further, the bucking of the travelling beam at a supercritical speed restricts 
its speed of operation. Thus the objectives of the vibration engineers have been 

(i) to increase the critical speed, and 

(ii) to reduce the level of vibration in the subcritical speed regime. 

As mentioned earlier, the critical speed can be increased by mounting the pulleys 
flexibly on the foundation [37]. Other methods like increasing the band tension by 
heating or applying a belt-tensioner system have also been used [75]. A new method 
of stabilizing the band at critical speed by introducing a parametric excitation has 
been suggested recently [76]. This method is analogous to stabilizing an inverted 
pendulum by means of parametric excitation. Difierent control strategies, both passive 
a.nd active, are adopted to reduce the unwanted oscillations of a travelling system 
during operation. Each controller has its own advantages and disadvantages. While 
the passive controllers are easy to implement, their inability to operate outside certain 
frequency band-widths limits the application [77]. On the other hand, an active 
controller, though difficult to implement, can be used for controlling systems subjected 
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to broad-band excitations [78]. The vibration of. a travelling continuous system is 
passively controlled by letting them pass through a number of hydrodynamic bearings 
[79, 80]. As the damping force is sensitive to the length rather than the location of 
a bearing, it can be placed at any convenient position within the span [81, 82, 83]. 
Active controllers have also been devised by measuring the system-responses with the 
help of suitably placed sensors and applying necessary excitation (based upon these 
measurements) through actuators placed at proper locations [84, 85]. The required 
amount of excitation depends upon various control strategies used. The simplest of 
such strategies is known as ‘modal control’ whereby a few of the normal modes are 
actively controlled [86]- [88]. However, the adverse ‘control’ and ‘observation’ spill- 
over effects of the unmodelled modes may make the controller ineffective [89, 90]. 
For a travelling beam, the observation spill-over effect, which is more dangerous than 
the control spill-over effect, may be reduced by using a comb-filter having frequencies 
that are adaptively changed with the axial speed [91]. Frequency-domain methods of 
active control have also been tried to control the vibration of both beams and strings 
[92, 93]. A novel and altogether different approach of controlling the vibtration of a 
travelling string by absorbing the energies of the associated travelling waves at the 
boundaries has been proposed [94]- [96]. This type of control is known as ‘wave- 
cancellation’ or ‘boundary’ control and has been applied for many axially stationary 
systems [97, 98]. Since the corrective measures are taken at the boundaries, the spill- 
over effects are circumvented altogether. However, if the sensor is placed at any point 
away from the actuators, detrimental effects may be resulted. This is due to the finite 
time taken by the waves associated with the transverse vibration to travel from the 
observer location (wherefrom the information of the oscillatory motion is collected) 
to the actuator location (where the corrective measure is taken). The problem of 
the non-collocated’ controller has been removed by introducing a suitable time delay 


10 



between the sensor and the actuator [99]. Other types of controllers are also suggested 
to achieve specific objectives. As for example, a combination of feedback and feed- 
forward controllers has been attempted to solve the chain-sprocket mismatch problem 
which is caused by the boundary oscillations [100] . 

1.2.4 Free and Forced Vibrations of Systems Having Variable 
Span 

So far the discussion has been restricted to fixed-span travelling members. However, 
during vibrations of a robot arm having a prismatic joint or of a flexible appendage in 
space crafts, the length of the span and hence the inertia changes continuously. These 
systems are usually modelled as a flexible cantilever beam of variable length and the 
end load, if any, is considered as a tip mass [101]. The response of a beam with 
either monotonically or periodically changing length is obtained by various methods, 
namely, Galerkin’s technique using time-dependent assumed modes [101, 102], finite 
element method with time-dependent elements [103], analytical method [104, 105] by 
transforming the time-dependent boundary conditions to the time-invariant ones, etc. 
The periodic response of a beam with monotonically increasing or decreasing length is 
not possible due to the non-conservative nature of the system. The amplitude of the end 
point of the beam, moving with either uniform or non-uniform speed, increases during 
retraction (i.e., when the length is shortened) and decreases during extrusion (i.e., when 
the span lengthens) and the frequencies of oscillation change significantly. Even the 
presence of a tip mass or the flexibility of the channel which controls the movements 
of the beam does not alter the stability characteristics [106]. It has been proved 
mathematically that these stability characteristics are ubiquitous for any system having 
continuously expanding or contracting domain [107]. It should be mentioned that the 
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reverse trend obtained in reference [102] is the result of an incorrect formulation of 
the equation of motion and is erroneous. For a uniformly extruding beam, the flutter 
instability has been reported if the speed is increased beyond a certain value [106]. 
Similarly, for a uniformly retreating beam, the stability is regained beyond the same 
speed. Parametric vibrations have been observed in a beam having periodic axial speed 
[104, 105, 102]. The vibrations of a beam during simultaneous axial movement and 
rotation about a transverse axis have also been studied [108]. It is apparent from 
the above discussion that a travelling continuous system having a single support can 
be unstable quite often. However, very few researchers have addressed the problem of 
controlling these instabilities. A variable order adaptive controller has been proposed 
for this purpose [109]. 

1.2.5 Effect of Non-Linesirities 

The above discussion points out the richness of literature devoted to various aspects 
of the linear vibrations of travelling continua. Nonetheless, the effects of the geomet- 
rical non-linear , terms become substantial as the axial speed of the travelling member 
approaches the critical value. It has been noticed long back that the non-linear terms 
significantly affect the characteristics of the hyperbolic equation of motion of a trav- 
elling string [110]. Further, different phenomena observed experimentally cannot be 
explained from the linear theory. For example, the shape of the response envelope of a 
boundary-excited string travelling at a subcritical speed is identical to the first mode. 
The amplitude of vibration, however, increases continuously until the critical speed is 
reached, when the amplitude decreases suddenly and the higher modes are generated 
[111]. Furthermore, if a point harmonic forcing is applied within the span, the string 
then beyond a certain axial speed, starts oscillating in a direction perpendicular to the 
plane of excitation. These two phenomena known, respectively, as ‘mode-jumping’ and 
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‘ballooning’ cannot be explained by the linear theory [111]- [113]. In order to explain 
these experimentally observed phenomena, a non-linear model of the travelling string, 
neglecting the longitudinal oscillation, has been suggested [111]. More accurate models 
of non-linear string and beam, simply-supported at both ends, have been subsequently 
developed [114]. Free vibrations of travelling continuous systems have been studied 
following these models. Non-linear dynamics of a travelling cable has also been studied 
[115, 116]. The effects of the non-linear terms on the natural frequencies of axially 
moving members are obtained by using Galerkin’s technique with the linear normal 
modes of both axially stationary and travelling systems as the basis functions. This is 
followed by either a harmonic balance method or a perturbation technique [114, 117]. 
These studies reveal the hardening nature of the non-linearities. In reference [117], 
the analysis has been extended to the supercritical speed regime. 

The modal coupling during the non-linear oscillation restricts the analysis of the 
forced response to the near-resonance conditions. Consequently, the near-resonance 
response of a travelling string, to a harmonic displacement excitation at the boundary 
with a frequency close to the first natural frequency, has been obtained both analytically 
and experimentally [118]. A perturbation method was used. The Galerkin’s technique 
followed by harmonic balance method has also been used towards this end [119, 120]. 
It should be pointed out that a one-term approximation in this method gives wrong 
answers since the effect of the gyroscopic term is not included. The ‘modal jumping’ or 
‘ballooning’ has not been considered in the above references. Although various authors 
deal with the latter phenomenon [112, 113], the explanation of the former has been 
attempted only numerically [111]. 

For a travelling beam having a single fixed support, different approaches are taken 
to derive the non-linear equation of motion. A modified version of the Hamilton’s prin- 
ciple, suitable for the systems with changing mass has been used by several researchers 
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[121]- [124]. A new way of formulating the equation of motion, within the framework 
of Kane’s dynamical theory, has also been reported [125, 126]. In this method the uni- 
form axial motion of the beam is taken as a constraint [127, 128]. The free vibrational 
analysis of the non-linear equation of motion shows that the non-linear terms do not 
change the stability characteristics predicted by a linear analysis [129]. The effects of 
the non-linearities are felt in the variation of frequency-contents of the time response. 

However, as mentioned in Section 1.1, the non-linear terms may change the dynamic 
behaviour drastically from that predicted by linear theories. The major qualitative 
changes can be summarized as follows: 

(i) The non-linear terms couple the vibrations of the two spans in a serpentine belt- 
drive which are assumed to be decoupled in the linear model [130]. 

(ii) The harmonic edge loading, meant to excite the torsional oscillation in a travel- 
ling beam, introduces large transverse vibration [131]. This phenomenon, known as 
‘autoparametric excitation’ ( or ‘internal resonance’) is possible only due to non-linear 
coupling between the transverse and torsional oscillations. 

(iii) Contrary to the prediction by a linear theory, the non-linear term restricts the 
exponential growth of the vibration amplitude of a parametrically excited string to 
‘limit cycles’ [132]. 

(iv) The analysis of free vibration neglecting the non-linearities leads to the prediction 
of flutter in a simply-supported beam travelling at a supercritical speed [16]. This 
prediction is, however, wrong [133]. Although the proof has been given for a simply- 
supported pipe, this is true for a travelling beam as well. 

(v) Chaotic response has also been observed for a parametrically excited travelling 
beam [76]. This phenomenon can never be explained by a linear theory [134], 

Prom the above discussion, it is apparent that the omission of the non-linear terms 
can lead to erroneous results. Although some progress has been made towards under- 
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standing the effects of non-linearities on the oscillatory motion of a travelling string, the 
study is far from being complete so far as the travelling beams are concerned. Finally, 
the controllers proposed for the travelling continuous systems should be re-examined 
to adjudge their efficacies in the presence of the non-linear terms. 

1.3 Objective and Scope of the Present Work 

The major objective of the present thesis is to inquire into the effects of the geometric 
non-linearities on the vibration of an Euler-Bernoulli beam moving over two simple 
supports. Various operating conditions generating different types of excitations are 
considered. The non-linear equation of motion is derived in Chapter 2. Since the 
small non-linear term in the equation of motion is considered as a perturbation to the 
linear equation, a detailed knowledge of the linear system is essential. The rest of 
Chapter 2 is devoted to study the free and forced linear vibrations by two apparently 
distinct methods, namely the ‘modal analysis’ and the ‘wave-propagation method’. 
While the discussion of the former is limited to the existing results, the latter method 
is developed for the first time. The numerical results obtained from these two methods 
are compared. Damping is neither considered here nor will be considerd later except 
during the occasions where the presence of a damping term changes the qualitative 
behaviour of the system. 

In Chapter 3, a simple and computationally efficient method of obtaining the near- 
resonance response of a non-linear continuous system is presented. The method is 
elucidated, for the sake of simplicity, with an axially stationary beam having either 
simply-supported or clamped-clamped boundary conditions. The non-linear normal 
modes of the beam, similar to the linear normal modes, are first derived using a per- 
turbation as well as the wave-propagation method. These non-linear normal modes 
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axe subsequently used to obtain the near-resonance response. The results so obtained 
are then compared with experimental results available in the literature. The efficacy 
of the present method over the Galerkin’s technique is highlighted. 

The results of Chapters 2 and 3 are extended in Chapter 4 to derive the response of 
the non-linear travelling beam under both resonant and hard non-resonant excitations. 
A concept of non-linear complex normal mode is first presented. The near-resonance 
response has been subsequently obtained using these non-linear complex normal modes. 
The methods of deriving both the non-linear complex modes and the near-resonance 
response resemble closely the procedure outlined in Chapter 3. Like in any other non- 
linear system multiple values are obtained for the near-resonance response. However, 
all of these values are not physically realisable. The physically observable results are de- 
lineated by carrying out a stability analysis. When the excitation is hard non-resonant 
one, the modal method encounters difficulties. The method of wave-propagation, de- 
rived for the linear system in Chapter 2, has been effectively extended for such a 
situation. 

Chapter 5 deals with the effects of the non-linear term on a parametrically ex- 
cited travelling beam. The source of the parametric excitation is taken to be the 
tension fluctuation due to the imbalance in one of the pulleys which is mounted on an 
axially-flexible support. Single or multiple limit cycles, generated by the combination 
of non-linearity and parametric excitation in the presence of small viscous damping, 
are obtained using the non-linear complex normal modes, The stable limit cycles are 
separated from the unstable ones by carrying out a stability analysis. The response 
of the beam to simultaneous parametric and external harmonic excitations is also ob- 
tained. The non-linear complex normal modes are once again used for that purpose. 

In this case too, a stability analysis to determine the physically observable solutions is 
presented. 
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The stability of an accelerating or decelerating non-linear beam is presented in 
Chapter 6. This situation arises during the starting or stopping phases. A Lyapunov’s 
stability analysis is first carried out to study the stability characteristics of the beam 
having any arbitrary acceleration. For a decelerating beam, the stability is determined 
using a perturbation technique based on Multiple-Time-Scale (MTS) method. The 
uniform deceleration is taken as a small perturbation to the uniform axial speed. 

A controller to reduce the vibration of a non-linear beam moving with uniform 
axial speed is discussed in Chapter 7. A passive controller in the form of a roller guide 
is proposed. The finite compliance of the guide in the transverse direction, and the 
frictional force appearing at the guide-beam interface and acting in the axial direction 
are incorporated. The results of the linear analysis, for both free and forced vibra- 
tions, are used to obtain the non-linear complex normal modes that are subsequently 
utilised to get the near-resonance response. The effectiveness of the roller in reducing 
the amplitudes of both linear and non-linear vibrations of the beam under harmonic 
excitation has been brought out with the help of numerical examples. 

The main conclusions of the present work are summarized in Chapter 8. Some 
directions for further research are also mentioned therein. 
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Chapter 2 


LINEAR VIBRATION OF A 
TRAVELLING BEAM 


2.1 Introduction 

First, the equation of motion for the transverse vibration of a simply-supported trav- 
elling beam including the non-linear term is derived in Section 2.2. Considering the 
physical nature of the problem, the equation of motion has been sindplified from the 
more accurate results reported by Thurman and Mote [114]. For a slender beam, the 
non-linear term can be taken as a small perturbation to the linear equation. There- 
fore, before embarking on the non-linear analysis, the linear anlysis is presented in this 
chapter. 

The linear free vibration of a travelling beam has been studied by various analytical 
techniques that are suitable for solving a differential equation with given boundary 
conditions [16], [37]- [42]. The major objective is to determine the natural frequencies 
and mode shapes. The connection between the vibration of and the wave propagation 
in continuous systems provides better understanding of the oscillatory motions. It is 
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well known that in an axially stationary string, two waves, travelling in the opposite 
directions at the same speed, are reflected at the boundaries to generate the normal 
modes of vibration. The free vibration of an axially stationary Euler-Bernoulli beam 
has also been studied using the wave-propagation approach. A ‘phase-closure principle’ 
has been derived [135] that shows that the normal modes correspond to a phase change 
of an integer multiple of 2-k for the propagating as well as the evanescent waves, as 
they return to their starting point after traversing to and fro once across the span. 
The phase-closure principle still remains valid during the modal vibrations of an axially 
moving string [9], where the speed of the wave propagation does not remain equal in the 
upstream and downstream directions. Owing to the dispersive nature of the medium, 
the wave propagation in an axially moving beam is more complicated. Consequently, 
only an approximate, rather than the exact, result has been obtained for the natural 
frequencies of such a beam [40]. 

Forced vibrations of a travelling continuous system are usually studied using Galerkin’s 
technique [58] or modal analysis [16]. In these methods, the equation of motion (which 
is a partial differential equation) is discretized into a set of countably infinite number 
of uncoupled ordinary differential equations. This is achieved by assuming the spatial 
distribution of the response in terms of various orthogonal (or ‘basis’) functions. In 
practical situations, only the first few terms corresponding to the lower modes are con- 
sidered. The error due to this modal truncation can be eliminated by constructing a 
closed-form transfer function (i.e., the Laplace transform of the Green’s function) by 
analytically solving the boundary value problem [60, 61]. An interpretation of the 
transfer function of an axially moving string, in terms of the waves propagating in two 
opposite directions [63] without distortion, simplifies the computation of the response 
to any displacement excitation. However, in a beam, where the wave propagation gets 
significantly complicated, the determination of the response by the wave-propagation 
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theory is difficult. 

The free and forced vibrations of a travelling beam are studied using modal analy- 
sis and the wave-propagation theory in Sections 2.3 and 2.4, respectively. The details 
of the modal analysis, easily available in the open literature, are omitted. But the 
method based on wave propagation is carried out in details. The frequency equation of 
the free modal vibration of the beam is derived using the phase-closure principle. Since 
the characters of the waves change drastically as the axial speed is increased beyond 
a certain value, the phase-closure principle is applied differently in these two speed 
regimes. The justifications of some approximate methods for obtaining the natural 
frequencies are also discussed. In the presence of an external forcing, the closed-form 
transfer function is first obtained using the wave-propagation theory, which is then 
used to determine the response to any arbitrary excitation. Numerical results obtained 
from the above two methods are compared by evaluating the response to a point har- 
monic excitation. The efficacy of the wave-propagation method for determining the 
forced response is clearly brought out. However, it should be mentioned that both the 
methods are important for carrying out the non-linear analysis presented later. The 
modal analysis is more suitable when the beam is resonantly excited whereas the wave- 

propagation method is advantageous if the beam is subjected to a hard, non-resonant 
excitation. 


2.2 Equation of Motion 

In this section, the equation of motion of the undamped, planar, transverse vibration 
of a slender beam, (see Figure 2.1) having an initial tension and moving axially 
with a uniform speed c* between two stationary and widely separated frictionless sup- 
ports, IS derived retaining the geometric non-linear terms. To this end, the following 
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Figure 2.1: Schematic diagram of a travelling beam. 


assumptions are made: 

(i) The amplitude of the longitudinal vibration, is order of magnitude smaller 

than the amplitude of transverse vibration, t). Mathematically, t) = 0(w*(^, 
where ^ represents the distance measured along the direction of axial movement, in a 
stationary reference frame, and t denotes the time. 

(ii) The radius of gyration of the beam cross-section, r, is much smaller compared to 
the span length I (i.e., r Z). 

(iii) The frictionless supports prevent the longitudinal oscillation and act as simple 
supports so far as the transverse vibration is concerned [114]. 

With these assumptions, the equations of motion of the coupled longitudinal and 
transverse vibrations of the beam are [114]: 
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where the area of the beam cross-section is A with second moment about the neutral 
axis Iz — Ar^, p and E are, respectively, the density and Young’s modulus of the beam 
material. 

The boundary conditions are given by 
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equations (2.1)-(2.5) can be rewritten, respectively, as 


d^u d^u 
dr^ ^dxdr 


+ c' 


dx^ 




7^ dx^ 


&^w d^w , o ^ ^ d'^w d^w o d 


dr^ dxdr 


^ ^ 72 ^ 

\j'^^dx dx 2^dx 


and 


u(0,r) = it(l,r) = 0, 
w {0, t ) = w (1, t ) = 0 


d^w(Q,T) ^w{ 1 ,t) 
dx^ dx^ ^ “■ 


( 2 . 6 ) 


(2.7) 

, ( 2 . 8 ) 
(2.9) 
( 2 . 10 ) 

(2.11) 
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The non-dimensional axial speed, c, has been obtained using the longitudinal wave 
speed and the slenderness ratio 7 . The significance of the inclusion of 7 will be apparent 
during the non-linear analysis in terms of the small parameter ( 7 ^/ 2 ). 

For a slender beam, i.e., with 7 -C 1, one can neglect the longitudinal inertia term 
in equation (2.7) and gets (neglecting terms of order higher than 7 ^) the following 
equation: 


1 d'^u _ 

7 ^ dx^ dx dx^ ’ 
which, when integrated twice with respect to x, yields 

u{x,T) = ~j^ dxi+xJ[{T) + f!i{T). (2.12) 

Ensuring u(x, r) satisfies the boundary conditions given by equation (2.9), the con- 


stants of integration /((r) and /^(t) are obtained as 

7“* 


/'(r)=Oand/;(r) = l-^(^)»di. 


Now substituting equations (2.12) and (2.13) into equation (2.8), one gets 


d'^w 


d'^w 


dr^ ^‘^^dxdr 


d^w d^.w 




= e 


Jo ^dx^ 


dx 


dx^ ’ 


(2.13) 


(2.14) 


dx^ dx^ 

2 

where e(= ^) is a small parameter, i.e., e -C 1. Finally, equation (2.14) together with 
the boundary conditions ( 2 . 10 ) and ( 2 . 11 ) describe the non-linear transverse vibration 
of the travelling beam. 

The equation of motion in the presence of a non-dimensional external force /(x, r) 
is 


d'^w 

dr^ 


+ 2c 


dxdr 


+ (c? - r„) 


dx^ 


+ 


d^w 

dx^ 



+ f{x,T), 


(2.15) 


where /(x, r) = f*{x,t)l/{EA'y^) with f*{x,t) as the transverse external force per unit 


length at the location x. 
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For the gyroscopic system under discussion, equations (2.14) and (2.15) are often 
recast in the following state-space form [16] : 

A^-pBVF = eN (2.16) 

OT 

and 

dW 

-I- BVF = cN -f f, (2.17) 

or- 

where 



with K ={(? — ro)5^/5a;^ + d^/dx^, G = 2cd/dx and / as the identity operator. 

Both forms of the equations of motion, (i.e., the partial differential equations and 
the state-space form) are useful depending on the method of subsequent analysis as 
will be seen later. 

In this chapter, the linear analysis of the free and forced vibrations (i.e., with e = 0 
in equations (2.16) and (2.17)) are presented. In what follows, the free vibration is first 
studied using both modal analysis and wave-propagation approach. The response of 
the beam to external force excitations is next obtained again using both these methods. 


2.3 Free Vibration 

As mentioned in Chapter 1, the transverse waves in a travelling beam move with differ- 
ent speeds in opposite directions. Consequently, no standing (stationary) mode exists 
even for the linear vibration. However, the situations still exist when all the points 
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of the beam undergo harmonic oscillation with same frequency. These frequencies are 
called natural frequencies of the travelling beam. These natural frequencies and the 
corresponding complex normal modes can be obtained by either direct modal analysis 
or wave-propagation method as detailed below. 


2.3.1 Modal Analysis 


The linear normal modes of a travelling beam differ from those of an axially stationary 
beam in the sense that in case of the former, the phase difference between the oscillatory 
displacement and velocity at any point depends on the axial coordinate x. This implies 
the importance of considering both the displacement and velocity at a point during 
the modal responses of a travelling beam. This fact is taken into account by treating 
the displacement and velocity as two independent quantities, which can be done by 
recasting the equation of motion in the state-space form (i.e., equation (2.16) with 
€ = 0) as shown below : 

A^-l-BW = 0. (2.18) 

or 

For a harmonic solution, the response vector W {x, r) can be assumed as 


W(x,r) = 


(2.19) 


where i = y/^ and $ = 



The bar at the top indicates the complex conju- 


gate. 

Substituting equation (2.19) in equation (2.18) and equating the coefficients of 
and separately, from both sides, one obtains the following equation and its 


complex conjugate ; 


= 0 , 


( 2 . 20 ) 



where the partial derivatives in A and B are replaced by total derivatives. It can be 
verified that equation (2.20) is equivalent to the following equation : 




and the boundary conditions are 




( 2 . 21 ) 


(2.22) 


A set of values of the natural frequencies, lo^ and the corresponding complex normal 
mode shapes, ^(a:), can be obtained by solving equations (2.21) and (2.22). Assuming 
in the following form 

^(^) = E 

i=i 

and substituting it in the boundary conditions the following frequency equation is 
finally obtained: 

1111 

kl kl kl 

^ik2 ^ik^ 

kle^^^ kle^^^ kle'^* 

where kjS (j = 1, 2, 3, 4) are the roots of the algebraic equation 




-{Jf - 2<xo^k - (c2 - To)e + k^ = 0. (2.24) 

The ra-th natural frequency and the corresponding normal mode, denoted by and 
4>n, respectively, are obtained by (numerically) solving equations (2.23) and (2.24) 

simultaneously. It can be shown that for the special case To = 0, equation (2.23) takes 
the following simple form 


<?cS^ cos c - cci 2 cosh |] + [c^ _ 8(„i)2] sin | sinh | = 0, 
^ J ^ z 


(2.25) 
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if c < 2^'^ and 




o S 2 

cos c — cos — cos — 
2 2 


+ 


4c2u;i - - 
” 2 


. (T . 62 _ 

sin — sin — = 0, 
2 2 


(2.26) 


if <c< (Ccr)i(= tt), where 5i = - c^, <52 = -4w^ and a = 

The linear analysis remains valid so long as the axial speed is less than the lowest 
critical speed, (ccr)i- For a tensioned beam [114] (Ccr)i = y/ToT^. 

Attention may be drawn to the fact, that the frequency equation for c < was 
also obtained in reference [39] (see equation (7) therein). However, that equation, in 
terms of the symbols used in the present work, turns out to be 


c^cr^i 


<T 5i 

cos c — cos — cosh — 
2 2 


O’ S 

— 16ti;^ sin — sinh -^ = 0. 


This equation is erroneous and should be replaced by equation (2.25). 

The well-known orthogonality relations [16], satisfied by the complex normal modes 
<^„’s, are given by the following complex inner products: 


and 


[ A$n da: = 0 for all m and n, 
Jo 

— T* 

/ da: = 0 for all m ^ n. 

Jo 


Using the expansion theorem [18], any non-dimensional vector V - 

[ ^2ix) J 

fying the boundary conditions can be expanded in terms of the complex normal modes 


ux{x) 


(2.27) 


(2.28) 


> satis- 


as 


V = 53(an^n + bn^n)i 


n=l 


where a„ and 6„ are obtained, with the help of the orthogonality relations, as 


^ So dx 

Otn — 1 T 5 

/o^$J^A$„dx 
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Figure 2.2: Variation of the first linear natural frequency with axial speed. To = 1.0. 
and 

_ Jo^^^AVdx 
” /o'¥jA$„da:’ 

2.3.2 Numerical Results and Discussion 

In this section, some numerical results obtained from the theoretical analysis presented 
in Section 2.3.1 are discussed. The variation of the first linear natural frequency {u{) 
with the travelling speed, c, is plotted in Figure 2.2. The zero value of at c = (Ccr)i 
implies the divergence or buckling instability. Figures 2.3 and 2.4 show the real and 
imaginary parts of the first and second complex normal mode shapes, i.e., ^i{x) and 
<j) 2 ix), respectively, for three values of the travelling speed. The real and imaginary 
parts are neither symmetric nor antisymmetric about the mid-point. 
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2.3.3 Wave-Propagation Analysis 

In this section, the free vibration of the travelling beam is re-examined using the wave- 
propagation theory. To this end, first the harmonic waves associated with the modal 
vibration of a simply-supported travelling beam are considered. The study is then 
extended to set up both the approximate and exact frequency equations. 


2.3.4 Wave Propagation and Reflection at a Boundary 


The equation of motion for linear oscillation is obtained by putting e = 0 in equation 


(2.14), as 


d^w d'^w 
dr^ ^dxdr 


+ {c^ - To) 


d^w 

dx^ 


-t- 


d^w 

dx^ 



(2.29) 


To study the harmonic wave propagation, the response of the beam is assumed as 


w{x, r) = (2.30) 

Substitution of equation (2.30) into equation (2.29) leads to same algebraic equation 
(2.24). Using Descartes’ rule [136], it can be seen that the roots of this quartic equation 
(in the wave number k) belong to one of the following categories : 

(i) one pair of real roots (with one negative (^i, say) and the other positive {k 2 , say)) 
with the other roots as complex conjugates {k^ and k^, say), 

(ii) four real roots with one positive (k 2 , say) and the other three negative {ki, and 
k^, say). 

In what follows, for the sake of algebraic simplicity, the analysis is carried on with 

Tq = 0. However, the nature of the wave propagation remains the same for all values 
of To. 
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For To = 0, the roots k/s {j = 1, 2, 3, 4) are obtained as 

ki = c /2 - + cV4, k 2 = c /2 + + 

= -c/2 + and A :4 = -c/2 - (2.31) 

Thus, the response is given by 

w{x, t) = (2.32) 


The wave numbers ki and ^2 correspond to the waves propagating along the down- 
stream and upstream directions, respectively, and will be denoted as Ai- and A 2 - waves. 
Their phase velocities CPi and CP 2 differ by 


CPi - CP2 = 



a;' 





— 



h 


k2 


— c. 


For c < 2ViJ, the other two waves i.e., the . 43 - and 4[4-waves can be called ‘evanescent 
waves’, since their group velocities become imaginary (implying no propagation of 
energy) as shown below : 




CG, = ^ = CGt = ^ = 2iy/w‘ - c^ji 


dk 


dki 


It should be noted that unlike in an axially stationary beam (c = 0) [135], in the 
present case, the phases of the evanescent waves change as they move along the beam. 
However, for c > 2\/a7, all the group velocities are real and both the Az- and 414 -waves 
become downstream propagating waves. The existence of such four propagating waves 
has been reported in reference [ 10 ]. Since the wave-propagation characteristics change 
drastically in the two velocity regimes mentioned above, the phase-closure principle is 
applied separately for these cases, namely for 0 < c < c^ and Cd < c < (ccr)i = tt with 
Cd — 2\/u?. But before applying the phase-closure principle, the reflection of waves at 
the boundaries with simple supports are briefly discussed. 
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Reflected wave 


I 



x = 0 


Figure 2.5: Reflection of wave at a simple support. 

Reflection of waves for 0 < c < c^: 

If the downstream propagating wave i.e., the .Ai-wave is the only wave impinging 
upon a boundary (see Figure 2.5), then the reflected waves will be both the upstream 
propagating {A 2 -) and evanescent (A 4 -) waves. The total displacement at a point is 

w{x, r.) = + Aae***® + ^e’***] e^'^ . 

which together with the boundary conditions 

, , d'^vj , , 

give the following relations 

Ax + A 2 + A 4 = 0 (2.33) 


and 

kfAi + 1^A2 “H k^A^ — 0. 

Using equations (2.31), (2.33) .and (2.34) one obtains 

^ (uj‘ — c/2^w< + c^/4) + ic/2xju}^ — (?/A ^ 

(w^ + cflxjiJ + <?IA) + ic/2Juj^ - c2/4_ 


(2.34) 


(2.35) 
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and 


X, = _ I . (2.36) 

_ (aj‘ + c/2^Ju}^ + c2/4) + icl2\fj^c^/ 4_ 

It can be seen that unlike in an axially stationary beam, the phase of the propagating 
wave does not change by tt after reflection at the boundary. The change of phase is 
given by 

cl2\j(J — c^ji 


Ciz = 7r+ l^tan 
= TT + Cl 


0(2^ b)^ — (P ‘1 A. 

UdlJ. — •' ' 

w' - c/2yJ(jj^ + c2/4 ■ + c/2‘^uj^ + 


— tan 


where 



- (cVif 
(u')2 + 2(c2/4)2 


(2.37) 


Reflection of waves for Cd < c<Tr: 

In this speed regime, the reflection can be studied by considering only the upstream 
propagating A2-wave impinging on the boundary. Three waves will be generated, 
namely the .Aj-, A 3 - and 244-waves. From the boundary conditions one gets 


(^41 -f 24$ + 244) = —A 2 . 

Since the 242-wave is the only wave reflected, the total phase change will be tt. The 
individual phase relationships between the 24i-, 243- and 244-waves need not be consid- 
ered here. 


2.3.5 Phase-Closure Principle 

In this section, the wave-propagation theory is used to obtain the natural frequencies. 
In what follows, for 0 <.c < cj, first the evanescent waves are neglected to obtain the 
natural frequencies approximately. Then the exact natural frequencies are obtained by 
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considering the evanescent waves as well. For ca < c < tt, there exists no evanescent 
wave and the exact natural frequencies are obtained. 

Natural frequencies neglecting evanescent waves for 0 < c < 

As explained in Section 2.3.4, the evanescent waves are characterized by their spatial 
exponential decay in the medium. Consequently, if the wave number of these waves is 
high, then the decay will be substantial and it is justified to neglect them altogether. In 
other words, the approximation becomes more accurate while determining the higher 
order frequencies. However, as can be seen from equation (2.31), the evanescent waves 
decay less with increasing axial speed. Therefore, the approximate derivation of the 
natural frequencies neglecting the evanescent waves is valid for low axial speeds. 

With the above limitations, the change of phase, as the propagating waves travel 
once around the span, is 

^^phase — 2(7r + Ci) + |A:i| + 1^121 , 

where the first term (within the parenthesis) in the right hand side signifies the change 
in the phase due to reflection at the boundaries and the other terras denote the phase 
changes as the propagating waves move along the beam. According to the phase-closure 
principle, for the n-th normal mode 

2(7r + €i) + + c2/4 = 2(n -f l)7r 

or 

€i -f- c2/4 = n-K. (2.38) 

Equation (2.38) together with equation (2.37) can be used to determine the natural 
frequencies. 

If, as a further approximation, ei is altogether neglected, the frequency equation 
becomes 

= n^TT^ - (2.39) 
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as obtained by Nelson [40], although from an altogether different approach. 
Calculation of the natural frequencies including the evanescent waves for 
0 < c < c<i 

In this section, the exact natural frequencies are derived considering the effects of the 
evanescent waves. For an axially stationary beam, the natural frequencies are obtained 
by considering the phase closure of either the propagating or the evanescent waves 
[135]. However, as already mentioned, in a travelling beam, the evanescent waves also 
contribute to the phase change during propagation. So, the phase-closure principle has 
to be used by considering both the propagating and the evanescent waves. 

Considering both forms of wave, the reflection at the boundaries can be analysed 
by considering the simple support at any point, say at a; = 0. The waves must satisfy 
the relations 

Ai -f- A 2 + A 3 -b A 4 = 0 (2.40) 

and 

kfAi + klA2 + klAz -b A:|A4 = 0. (2.41) 


Using these relations, one can see that if the Ai- and Aa-waves impinge on a simply- 
supported boundary, the resulting A 2 - and A 4 -waves are given by 


Similarly, 





Ai 

A 3 




1 1 

kl kl 

1 1 

u2 1,2 
^3 


-1 r- 


-1 r 


1 1 
t,2 

1 1 

kl kj 


Ai 

A 3 


(2.42) 



(2.43) 


Now, to use the phase-closure principle one needs to follow the waves until they reach 
their starting point after travelling to and fro once along the beam. This is carried out 
as explained below: 
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(i) Let Ai- and Aa-waves of respective strengths ai and as start from a; = 0 and travel 
towards a: = 1. If the strengths of the waves at a; = 1 become a'l and 03, respectively, 
one can verify that 


a, = e 


a '3 = e‘*^a3, 


or in the matrix form 


gifci 0 

0 e’*® 




, (say). 


(ii) The Ai- and Aa-waves reflect at the simply-supported boundary placed at a: = 1 to 
3deld the A2- and A4-waves of strengths 03 and 04, respectively. Since the phase change 
due to the reflection is independent of the location of the boundary, the transformation 
matrix is always given by equation ( 2 . 42 ). Thus, 


[BJ = - 

ju2 l2 ]L2 u 2 

^2 ^4 ^1 ^3 

. (iii) The upstream propagating (i.e., the A2-wave) and evanescent (i.e., the A4-wave) 
waves reach from a; = 1 to a; = 0 with strengths c 4 and a'^, respectively. The transfor- 
mation matrix is, 

0 1 r 0 1 

[^*3]= .1 = 


g-ifc2 

I 

0 



1 

0 

0 

1 

1 


0 

1 

— 1 

‘Iii 


(iv) The A2- and A4-waves get reflected at s = 0 and become the Ai- and Aa-waves 
with strengths % and a^, respectively. The transformation matrix is given by equation 
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(2.43). Thus, 


1 1 

kl kl _ 

The condition for the beam to vibrate in one of its normal modes can be stated as 

a'l = oi (2.44) 


m = - 


1 1 
kf kl 


and 

O3 = 03. (2.45) 

In terms of the transformation matrices, the phase-closure principle thus yields 




ai 


0-3 


fll 


= < 


O3 




(2.46) 


For non- trivial solutions of and 03, the following condition must be satisfied : 


det [[i?4][i23][i22][i2i]-[/]] = 0. 


(2.47) 


As shown in Appendix A, the above equation is same as the frequency equation (2.23). 
Calculation of the natural frequencies for Cd<c<7c 

In this speed regime, three waves, namely the Ai-, A3- and A.4-waves propagate in 
the downstream direction and the remaining A2-wave moves in the opposite direction. 
When a simply-supported boundary is placed at a point, say a; = 0, the four waves 
satisfy the following relations: 


Ai -1- A2 -f- A3 + A4. = 0 


and 


^A\ "h klA2 4" ^ 3^.3 k^A^ = 0 . 


If an A2-wave hits upon the boundary, the resulting Ai-, A3- and A4-waves are obtained 
from the above two equations. But since the number of unknowns exceeds that of the 
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equations, the strengths of the waves can not be calculated individually. However, their 
relative strengths AijA^ and A^fA^ are obtained as 

A1M4 1 11 ^ ^ I ^2M4 

Az/A^ J k\ k\ k^ k\ 1 

On the other hand, if the propagating Ai-, Az- and i44-waves get reflected from the 
boundary, the resulting >l2-wave can be obtained from the following relation: 

A2 — —(.^1 + Az + .A4). ( 2 . 49 ) 

It is to be pointed out that this motion is over-constrained. The waves must also 
satisfy the following equation due to the disappearance of the bending moments at the 
simply-supported boundaries: 

Aik^ + Azk^ "f" Azk^ -t- A^k"^ = 0 . ( 2 . 50 ) 

To get the normal modes, the propagation of the waves is to be followed as detailed 
below : 

(i) Let an ^la-wave, of strength 02 be reflected from the simply-supported boundary 

placed at a; = 0 . Three waves, namely, the Ai-, Az- and i44-waves, are generated. The 

strengths of these waves can not be determined uniquely, but the relative strengths 

01/04 and 03/04 can be obtained in the following matrix form : 

N r' n~^r "] f ^ ( 

O1/O4 11 111 O2/04 0,2! 0 >\ 

> = - < 

03/04 J [ ^1 ^3 J L *2 *4 J i 1 J I 1 

(ii) The above Ai-, Az- and yl4-waves of respective strengths Oi, 03 and 04 start trav- 
elling from a: = 0 and reach x = 1 with strengths o^, O3 and 04, respectively. It is seen 
that 

Oj = O 3 = e**^a3, = e‘*^‘‘o 4 . 
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Since only the values of the relative strengths 01/04 and 03 / 04 are known at x — 0, the 
corresponding values at x = 1 are obtained as 


I a'l/a; 
i 03/04 , 




0 


0 



01/04 

1 

. = [iJ,] . 

di/ d 4 

< 

1 ^ 


^ 03/04 J 


[ d 3 /d 4 ^ 


(say). 


(iii) The propagating Ai-, Az- and .r 44 -waves are reflected from x = 1 and produce 
an ^ 2 -wave of strength Oj. As equations (2.49) and (2.50) are to be satisfied at the 
boundary, one gets the following relation in terms of the relative strengths: 



> = 


n-l r 


k\ 


1 1 

k\ kl 


oi /04 

i J 


= [^ 3 ] 



(iv) Finally the Aa-wave of strength o^ reaches the point x = 0 with strength Oj. Thus, 


(I 2 




from which the following transformation matrix is constructed: 


02 / 04 1 

Q-iik2~k4) Q 


03/04 ' 

1 ■ 1 

► = [Ra] < 

O2/04 

1 r = 


i 


i 1 J 

1 

0 

t 


1 1 J 


[ 1 J 


According to the phase-closure principle 

ft 

0*2 ^ ^2? 


or 


J ^4 ““ Ot2 j d4 • 

In terms of the transformation matrices, the above requirement is equivalent to 


[[i24][i23][i?2][^ll]-/] 


1 


> = 0 . 
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For non-trivial solutions 


det [[i 24 ][i 23 ][ii: 2 ][i 2 i] - [/]] = 0. (2.51) 

One can easily verify, following the results given in Appendix A, that equation (2.51) 
is again identical to the frequency equation (2.23). 


2.4 Forced Vibration 

In this section, the forced vibration of the travelling beam is studied using the modal 
analysis as well as the wave-propagation method. 


2.4.1 Modal Analysis 

The equation of motion of a linear travelling beam, excited by an external force f{x, r), 
can be written, after substituting e = 0 in equation (2.17), as : 


+ BW = f. 

OT 


(2.52) 


The response W{x^t) and the initial conditions are expanded in terms of the complex 
normal modes in the form 


W(l,r) = i:(C„(r)«„+C„(T)4'„) 

n=l 


and 


00 


w'(^.o) = E(Cn(o)«.+<.(o)«„). 

n=l 


■ (2.53) 


(2.54) 


Combining equations (2.52) and (2.53) and the orthogonality relations given by equa- 
tions (2.27) and (2.28), one can easily get the following equation and its complex 
conjugate 

dCn /o'Cfd^ 

■57 - “»C» - 71^. 
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The above equation can be solved as an initial value problem after obtaining the initial 

conditions from equation (2.54) and using the orthogonality relations. In what follows, 

two kinds of forcing, namely, a point-impulse and a harmonic loading, are considered. 

For a unit-impulse excitation at point a; = Xq, i.e., f{x,T) = 6 {t—to)5{x—Xo) where 

d(a;) denotes the Dirac d-function, the values of Cn(7’) ’s with zero initial conditions can 

be obtained as r , _ i ■ 

^ , . [lo - ro)dri 

yielding the response as 


fo Mn ~ 2ic /o dx 




00 


W 


{x,r) = Yl 


(j>n{x)(pn{xo) 


n=i2uji^S^<f)n(f>ndx-2icS^<j)^^dx 


to) (2.55) 


and zero for r < tq, where c.c. denotes the complex conjugate of the preceeding ex- 
pression. 

The steady-state response (neglecting the effects of the transients) of the beam to. 
a harmonic loading f(x, r) = f{x) cos Qr can be written as 


W{x,t) = i [P{x)e^^ -h P(x)e-'"’‘) , 


(2.56) 


where 0 is the non-dimensional frequency obtained from the frequency of excitation, 
0.*, through the relation 



Following the modal analysis detailed above, one finally gets 


where 


S(Pn^n + 9n^n), 

fl=l 


(2.57) 
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and 


(2.58) 


Qn 




with fi = s 


i(n + Oj},) fo A$n dx ’ 

/(x)/2 ' 

0 

It is to be noted, that for the stationary beam <^n = 4>n fo = 

2(a;^)^/o ^^dx. Thus, the response expressed by equation (2.56), after expanding 
and fi, becomes 

W(x, r)=(f] \ cos Qt, 

which is well established in the literature [137]. 


2.4.2 Wave-Propagation Analysis 

The response of the travelling beam first to a point impulse- excitation and then to 
any arbitrary excitation are now derived using the principle of wave propagation. To 
this end, the following observations of the free vibration, presented in Section 2.3.4 are 
important. 

The natures of the waves generated during the free vibration are different in two 
axial speed regimes, namely, c < Cd and c > Cd- In the first case, one propagating 
(Ai-wave) and one evanescent (Aa-wave) waves travel in the downstream direction. 
Similarly, the upstream waves comprise of one propagating (4,2-wave) and one evanes- 
cent (4.4-wave) waves. In the second case, in place of the two evanescent waves, two 
more downstream propagating waves (43- and 44-waves) are generated. In either case, 
equation (2.23) is satisfied during a modal vibration. 

The easiest way to understand the forced response of a continuous system through 
wave-propagation approach is to analyse th,e equation of motion in the frequency do- 
main. This is done by writing the equation of motion in terms of the Laplace transform 
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of the variables, taken with respect to time, as follows : 


s^w+2sc^+{c^-To)^+^ = f{x,s)-s^{x,0)-w{x,0)-2c^{x,0), (2.59) 

where w{x, s) and f(x, s) are the Laplace transforms of w(x, r) and f{x, r), respectively. 
The initial configuration and velocity of the beam are. assumed to be trivially zero, i.e., 
at r = 0, w = 0, dwjdr = 0 and dwjdx = 0 for all x. Further, with f{x,T) = 
6{x — xo)5{r — To), the equation of motion becomes 


dw 


d'^w d^w 


H- 2.C- + (C^ - - x„)a-, 


(2.60) 


since f{x, s) = S(x - Xo)e~^'^. 

The particular integral of equation (2.60) can be obtained as 

j=l \m=l, mjij ^m) J •'0 

For x> xo, the particular integral becomes 


(2.61) 


ty(ar, s) = 

j=i 


where 


n 


, Q-sro 


A 1 ffc* - fc» r ’ 

1=1,771547 

and k^s are the roots of the polynomial 


(2.62) 


s'^ + 2isck*-{(^-To)ik*)^ + ik*)^ = 0. . (2.63) 


For any complex value of s, the roots of the above polynomial are either real or com- 
plex. Among the complex roots, those having positive -(or negative) imaginary parts 
correspond to the downstream (or upstream) evanescent waves. Of the real roots, the 
negative ones correspond to the downstream propagating waves and the positive ones 
to the upstream propagating waves. For s = icj, the possibilities of the two kinds of 
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waves are already mentioned in Section 2.3.4. Depending upon the relation between c 
and u, either two propagating or four propagating waves may exist. 

Prom the physical nature of the problem, it can be argued that all the terms in 
the right hand side of equation (2.62) are not present at a point x > Xq. At the point 
X = Xo, there exist four waves of strength Dj (j = 1,2, 3 , 4), out of which only those 
propagating in the downstream direction have an effect on a point lying in the same 
side of Xo- The effects of the upstream-going waves will not be felt at these points. For 
example, for s = iw (a> > 0 and c < ca) only the contributions of Ay- and Aa-waves are 
to be retained. But for c > c^, all the downstream propagating waves (i.e., the Di, Dz 
and £>4 terms) contribute to the response at a point a: > a;o. 

Thus, the omission or inclusion of any term of the right hand side of equation (2.62) 
depends upon the values of kys. In the following, the response is calculated only for 
s = ioj and c < Cd- The other cases may be studied in an analogous manner. The 
particular integral (see equation (2.62)) now becomes, 

w{x,s)'— ^ fora;> a:o, (2.64) 

i=i,3 

The particular integral for x < xq can be found out by first substituting x' = —x in 
equation (2.60). Following the same steps carried out for x > xq, and noting that —kys 
are the roots of the algebraic equation analogous to equation (2.63), the particular 
integral is finally obtained as 

4 

t&(x, s) = — for x < Xq. 

j=i 

Following the same reasoning, used for x > xo, the particular integral in this case has 
terms corresponding to the waves going only in the upstream direction. For the case 
considered here (s = iu and c< Cd), 

tD(x, s) = — 51 for x < xq . (2.65) 

J=2,4 
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It is to be noted that the particular integral gives the response of the beam only due 
to the direct influence of the force. However, the waves generated by the external 
excitation are reflected at the boundaries and new terms appear. These terms are 
taken care of by retaining the complementary function. The total response is thus 
given by 

w(s,s) = for x > xq, (2.66) 

i=i j=i,3 

= for x < Xq. (2.67) 

j=l j=2,4 

One can also verify that the continuity of displacement, slope and moment exist at 
X = Xq. The unknown constants C/s can be obtained by ensuring that the response 
satisfies the boundary conditions. This yields 

( 2 . 68 ) 

1 

) 

ik* 
e 4 

•xo) 0 

[ 0 (A:|)2e-**3(i-®o) q 

{^} = {C'lj C' 2 » C's, and {D} = {Di, —D 2 , Ds, The values of Cj’s become 

{C} = -[A-]-‘[Fl{£)} = -[A‘)^[i?){P}/det[A-J, (2.69) 


0 = [A'KC} + [F]{D} 


where 


[A*] = 


1 


ik* 

e 1 


[F] = 


1 


ik’^ 

e**2 


1 

(K? 

ik’^ 


(/:i)^e**i (A:2)^e**'* (AiD^e**^ 


0 

0 

xo) 


g— tfejXo 

(fe2)^e“**2®° 


0 

0 

-iktd- 
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where [A*]adj is the adjoint matrix of [A*]. Substitution of the values of C/s into 
equations (2.66) and (2.67), gives the transfer function for a travelling beam and is 
denoted by G{x, xq, s). 

The response of the beam to any arbitrary excitation f{x, t) is obtained from the 
Laplace transform of the response, given by 

w{x, s) = / G{x, C, 5)/(C, s) dC- 
Jo 

To retrieve the time-response w{x,t), to a unit impulse excitation, one has to take the 
inverse Laplace transform of w(x, s) as 

1 r^+iY 

w(x, t) = lim / w(x, s)e^'^ds, for r > 0 (2.70) . 

^ ' Y-^oo 27rtJ l3-iY ^ ^ ^ 

where the value of is so chosen that the integration converges. It is shown below that 
the integration can also be evaluated using the contour integration theory of complex 
variable. As seen from equations (2.62)-(2.67), the integrand contains a term for 
which the inverse transformation is carried out separately for r <tq and t > tq. ■ 

For T < To, the value of /3 is taken to be negative and the contour is devoid of any 
singularity (see Figure 2.6). Thus, 

w{x,t) = 0 for r < Tq. (2.71) 

For r > To, one must choose ^ to be an arbitrary positive number. Since the con- 
tour (see Figure 2.7) now contains countably infinite number of singular points, the 
integration can be performed with the help of the following two theorems [138, 139]. 

Theorem 1 : Let /(s) be a function which is analytic inside a simple closed path C 
and on C, except for finitely many singular points ai, a 2 , ■ • • Om inside C, then 

/ /(5)ds = 27ri53Res/(s)l , 

j=i 
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Figure 2.6: Contour of the integration (equation (2.70)) for r < tq. 
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Figure 2.7: Contour of integration (equation (2.70) for r > tq. 
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the integral being taken in the counterclockwise sense around the path C [138]. 

Theorem 2 : If /i{s) and / 2 (s) are analytic in the neighbourhood of Sq and if 
/i(so) # 0 but / 2 (s) has a simple pole at Sq, then the residue of /i(s)// 2 (s) at Sq is 
equal to /i(so)// 2 (so) [139]. 

The integration appearing in equation (2.70), when carried out, yields 


= S I a(diilw 1 




g“^n(T ro) _|_ g g ^ £qj. 2; > 

(2.72) 


where the value of the denominator is 


|(det[A-l)=|:Adet[Alf. 

The response cein be obtained for x < xq, after replacing Di and D 3 in equation (2.72) 
by — jDa and — U 4 , respectively. 

Before considering the steady state response of a travelling beam, it may be worth- 
while to mention that the response of a travelling string obtained by the wave-propagation 
method can be explicitly shown to be identical to equation (2.55). The details of this 
derivation is given in Appendix B. 

The computation of the steady-state response of a travelling beam to a harmonic 
point load is very convenient by the wave-propagation method. For a linear beam, with 
a harmonic excitation of frequency (i.e., /(x,t) = Fq5{x — Xq) cos Qr), the response 
is obtained as 

w(x,T) = ^G(x,xo,in)e^^^ + c.c.. (2.73) 

It may be noted that the response becomes infinitely large if the beam is resonantly 
excited with f2 = as det[A*] = 0. 
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2.4.3 Numerical Results and Discussion 


In this section, numerical results of the linear response, obtained by both the wave- 
propagation method and modal analysis are presented. The initial tension and the axial 
speed of the beam is taken as To = 1 and c = c/(ccr)i = 0.5 where (ccr)i = \/To -t- tt^. 

Figure 2.8 shows the variation of \G{x^ with the frequency of excitation 

(fi), when an exciting force of unit amplitude is applied at xq = 0.3 and the response 
is measured at a; = 0.5. As shown in the figure, the response amplitude is infinitely 
large if the excitation frequency, iQ, coincides with one of the natural frequencies of 
the beam and the point of either excitation or observation does not fall on the nodal 
points of the corresponding mode. Thus, the response amplitude at as depicted 

in the figure, is small since the point of measurement coincides with the nodal point 
where the response envelope corresponding to the second mode attains its minimum 
value. A large response is obtained for this excitation (i.e., with Q, ~ at any other 
point. Figure 2.9 shows a large response obtained at a: = 0.75. It must be pointed 
out, that an interchange between the positions of excitation and observation leaves the 
response unaltered. This reciprocity relation has indeed been observed in the numerical 
computaion. 

The response can also be computed in terms of the normal modes, whereby a 
solution is constructed in a series form. However, only a few terms of the infinite 
series are usually taken into account. Figure 2.10 shows the displacement profiles, 
obtained by considering only the first two linear modes of the beam, with an external 
force Fq = 13.22 and Q. = 0.9a;[ applied at Xq = 1/3. The four curves in the figure 
correspond to different time instants during one cycle of excitation. The continuous 
phase difference along the length of the beam is obvious in this figure. 
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Figure 2.8: Linear frequency response function to a point excitation at xq = 0.3 mea- 
sured at a: = 0.5. d = 0.5. 


asimAL tmm 

1.1. T., rAWFUl 


51 






(2 ‘r).-H 


53 


5 beam at four i 




For the excitation frequency close to a natural frequency of the beam, the corre- 
sponding normal mode contributes most significantly. But if the excitation frequency 
is away from any u;|^, the effects of all the neighbouring modes become significant. For 
a high frequency of excitation a large number of modes are to be considered, thus mak- 
ing the computation quite expensive. However, in the wave-propagation method, the 
response is obtained accurately in the closed form. Figure 2.11 shows the comparative 
results obtained by the wave-propagation method and also by the one- and two-term 
modal approximations. For the chosen frequency range, one-term approximation is 
very crude while addition of another mode makes the result closer to the exact. It 
has been found that for a frequency even the first mode approximation yields 

quite accurate results. 



Figure 2.11: Linear frequency response function to a point excitation at Xq = 0.3 

measured at a: = 0.5. cf = 0.75. ; wave propagation analysis; : one-term 

modal analysis; : two-term modal analysis. 



Chapter 3 


NON-LINEAR NORMAL MODES 
AND NEAR-RESONANCE 
RESPONSE 

3.1 Introduction 

The linear normal modes of a travelleing beam have been derived in Chapter 2. These 
modes have been subsequently used to obtain the forced response. In the non-linear 
case (i.e., e 0), such a modal vibration is not possible due to the coupling of the 
linear modes. Consequently, the analysis of the forced vibration gets complicated. This 
difficulty is nothing special for the vibration analysis of a travelling member, but ap- 
pears in almost every non-linear, multi-degrees-of-freedom system. However, periodic 
solutions exist for the equation of free vibration of even such a non-linear, multi- 
degrees-of-freedom system. During such a motion all the co-ordinates pass through 
their equilibrium positions simultaneously and also attain the maximum values simul- 
taneously. This fact, which is true during the modal vibration of a linear system, has 
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been used to introduce the concept of a ‘non-linear normal mode’, and a ‘non-linear 
natural frequency’. But generally the ‘non-linear normal modes’, unlike the linear ones, 
are not orthogonal to each other. Further, the shapes of the non-linear normal modes 
as well as the associated natural frequencies are not system properties. Both of these 
depend on the extent of motion. 

The non-linear normal modes of a raulti-degree-of-freedom system have been ob- 
tained originally using geometric methods [140]. The efficacy of the simple Harmonic 
Balance Method followed by numerical analysis for obtaining the non-linear normal 
modes of a multi-degree-of-freedom system has been brought out by Stupnicka [141] 
with the help of numerous examples. Numerical computation of near-resonance re- 
sponse of such systems, using the non-linear normal modes, is also amply demonstrated 
in the last reference. The concept of non-linear normal modes for an axially stationary 
beam (infinite degrees of freedom) with linear boundary conditions have been elabo- 
rated by several researchers, like Shaw and Pierre [142], Nayfeh and Nayfeh [143], 
King and Vakakis [144]. Different methods, namely, the invariant manifold theory 
[142], the energy method [144] and the multiple-time-scale method [143] have been 
used to achieve the same goal. 

In Section 3.2, a simple method of obtaining the non-linear normal modes of a 
continuous system is demonstrated with the help of an axially stationary beam, either 
simply-supported or clamped at both ends and having negligible longitudinal inertia. 
Towards this end, a combination of harmonic balance and a perturbation technique, 
instead of using a numerical approach following harmonic balance [141], has been 
used. Thus, the effects of non-linearities can be very easily tracked to different orders 
of approximation. The non-linear normal modes are subsequently used in Section 
3.3 to obtain the steady-state, near-resonance response for both symmetrically and 
asymmetrically located harmonic loads. The results obtained by the present method 
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are compared with those reported in reference [146]. In reference [146], the Galerkin’s 
technique was used for response computation and the results were compared with those 
obtained experimentally. 

In Sections 3.4 and 3.5, the non-linear normal modes and the near-resonance re- 
sponse are recalculated, respectively, taking the longitudinal inertia into consideration. 
The negligible effect of the longitudinal inertia on the transverse vibration of the beam 
is clearly brought out. 

Thus, the major objectives in this chapter are: 

(i) To develop a method for studying the free and forced vibrations of a non-linear 
continuous system. 

(ii) To provide a justification for neglecting the longitudinal inertia term while deriving 
the non-linear equation of motion of the travelling beam (see Section 2.2). 

Although the results are obtained only for an axially stationary beam, the method- 
odology is quite general and will be used again in Chapter 4 to analyse the free and 
forced vibrations of a non-linear travelling beam. 


3.2 Non-linear Normal Modes Neglecting Longitu- 
dinal Inertia 


In this section, the non-linear normal modes of an axially stationary beam without 
initial tension are obtained neglecting the effect of longitudinal inertia. The non-linear 
equation of motion of such a beam, obtained by substituting c = 0 and To = 0 in 


equation (2.14), is 



(3.1) 
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The n-th linear normal mode and the corresponding natural frequency (n = 
1,2, 3, • • •) are obtained from the above equation after substituting c = 0. To obtain 
the non-linear normal modes, the response w{x,t) is assumed as 


w{x, r) = cos Wnr. 


(3.2) 


Substituting equation (3.2) into equation (3.1) and equating the coefficients of cos tOnT 
from both sides, one gets 




d Vn 


j A = Tca 

dx^ 4 


Jq^ dx^ 


dx 


dV„ 

dx2 • 


(3.3) 


Now, for small non-linearity (i.e., e <C 1) the n-th non-linear mode "ipn and the associ- 
ated non-linear natural frequency cUn are assumed, respectively, as 


and 


‘lpn=4>n + ^ + ’ • ' , 


= (Wn)^ + + 


(3.4) 


(3.5) 


Substituting these expressions into equation (3.3) and subsequently equating the coef- 
ficients of the like powers of e from both sides of the resulting equation, one obtains 


-(^n)Vn + ^ = 0, 


dVn 


dx^ 


Jo ^dx^ 


dx 




(3.6) 


. (3.7) 


Equation (3.6) is trivially satisfied. To solve equation (3.7), (5^^^ is written in the 
following series form : 

«>= t 4‘Vi. (3.8) 




since the solvability condition [145] of equation (3.7) demands that cannot have 
any contribution from ij>n- 
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(3.9) 


Further, one can rewrite the right hand side of equation (3.8) as 


3 °° 




J=1 


where 


= 
3 


(•^0 (^) (-^o 

(/o ^?da:) 


(3.10) 


First equations (3.8)-(3.10) are substituted into equation (3.7), then equating the co- 
efficients of ^jS from both sides of the equation, one gets 




(3.11) 


and 


xn - 2„2 v' 


(3.12) 


Thus, the non-linear normal modes and the associated frequencies are obtained as 

(3.13) 


and 


I j* I ^ 2 ^3 

i>r, = 4>^ + W E r, , ,„T 

^ j^r> [(l.>i)“ - («y^J 




(3.14) 


The amplitude dependence of the non-linear mode shapes and the natural frequencies 
are clearly exhibited by the above two equations. 

Similar to the linear case, the phase-closure principle holds good even during the 
non-linear modal vibration. This is shown below by recalculating the non-linear normal 
modes and the natural frequencies using the wave-propagation theory. 


3.2.1 Wave-Propagation Approach 

To derive the non-linear normal modes, i.e., the shape of the beam which vibrates at 
a single frequency w, the response w(3:,t) given by equation (3.2), is substituted into 
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equation (3.1) to get 



after neglecting the higher harmonics. For a normal mode to exist, tj){x) is such that 
the phase of any wave should change by an integral multiple of 27r by travelling once 
to and fro along the beam. In the absence of non-linearity (i.e., with e = 0), the 
phase closure occurs when u> = and ip = (pn. In general, in the presence of the 
non-linear term, different linear modes are simultaneously excited. Hence, the normal 
mode corresponds to the frequency at which the phases of the waves corresponding to 
all the linear modes are closed simultaneously. Thus, to get the normal oscillation, one 
has to search for such a frequency. To this end, the following observation from the 
linear analysis is important. 

For the free vibration at any linear mode, say the n-th mode, one can write 

(Pr^ix) = + Aae**'* + + ^46^*, 


where the wave number k = \ju^ and the coefficients .41,^2, Az and A 4 bear a constant 
ratio amongst each other. Thus, four (two propagating and two evanescent) waves of 
definite wave numbers travel once across the beam in such a way that the phase of 
each wave gets closed after a time interval 27r/a;Jj. However, for a frequency different 
from ujI, say Cl, the same waves can be closed as they traverse once across the beam 
only by applying a suitable external force given by 


/» = [(“»)^ - 


(3.16) 


With the above forcing, the wave numbers and the relative strengths of the constituent 
waves do not change, but the phase velocity of each wave is changed to 


CP = - = -r--T = CPx-^ 
k k (jj\. Cl 
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where CP is the phase velocity of the wave propagating in the n-th normal mode and 
having the wave number k. 

Now returning to the non-linear case, for weak non-linearity one can write 




(3.17) 


and 


m=l 


m j 


(3.18) 


where o„ = a and = 0(e) for m^n. 

It can be said that the force arising out of the non-linear term (i.e. the right hand 
side of equation (3.15)) is just suitable for simultaneous application of the phase closure 
for all the waves participating in various linear modes. Keeping in view equation (3.16) 
for a linear normal mode and replacing <f>n{x) by aip(x) and by a; one can write 

o-mK = (-^f dx (3.19) 

Applying equations (3.17)-(3.19) together with the orthogonality relationships among 
the linear modes one gets the following results neglecting terms O(e^): 

(/o‘^«^'i^)(j'o'(^)'d-) 

!o<f>ldx 


= {^hf - 


(3.20) 


and 


3 3 

o„i = -a c 




The above equations yield results identical to those given by equations (3.13) and 
(3.14). 

In what follows, the results for simply-supported and clamped-clamped boundary 
conditions are explicitly obtained. 
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3.2.2 Simply-Supported Beam 


For the simply-supported non-linear beam, the boundary conditions are given by equa- 
tions (2.10) and (2.11), which are reproduced below: 


w{0,t) = iy(l,r) = 0, 


d^w{0,T) _ 5^w(l,r) 
dx^ dx"^ 


In this case, 


<pn = sin mrx and (w^)^ = (utt)^. 


(3.22) 


Substituting equation (3.22) into equations (3.10), (3.13) and (3.14), one gets 


=0 if i # n and = —{mr)^/2 


and 


= (nTT)^ 


, 3 2 

1 + -ea 


^„(a;) = (f>n{x). 


(3.23) 


Thus, for a simply-supported beam, the non-lineaf normal modes are the same as 
the linear normal modes and hence also orthogonal to each other. The above results 
are well established in the literature [143]. 


3.2.3 Clamped-Clamped Beam 


In this case, the boundary conditions are 

/« \ f-, \ n dw(0,T) dw(l,T) „ 

»(0,r)=Ml,T)=0, ^=0. 

The linear modes and the natural frequencies are given by 

(sin Hn — sinh //„) 


— 


(sin HnX — sinh finx) — 


(cos fin — cosh fin) 


(cos flnX — cosh flnx) 


(3.24) 
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where /^n’s are the roots of the transcendental equation cos cosh Hn 
the following relation is satisfied between the normal modes. 

P' dVn 


r d‘0„ 

Jo ^ ° if J + n = odd. 


= 1. Further, 


(3.25) 


Combining equations (3.13), (3.14), (3.24) and (3.25), one finally gets 


^ j+n=even,i?4n [(wj)^ - (n;^)2] 


(3.26) 


where aj^^’s are obtained from equation (3.10). 

From equation (3.26) one can see that the odd order non-linear normal modes 
comprise of only odd order linear normal modes and even order non-linear normal 
modes are given by a combination of only even order linear normal modes. Hence the 
odd and even order non-linear normal modes are orthogonal to each other. 

The non-linear normal modes for the clamped-clamped beam and the corresponding 

frequencies, taking only the first three linear normal modes into consideration, are 
obtained up to o(e) as 


V-i =<^i+0.006675€aV3, 

[l + 0.2351ea2] , 

(3.27) 

V '2 = <f> 2 , cjj = (4)2 [l - 1 - 0.41672ea2] , 

(3.28) 

^3 = <f>o- 0.05017ea2^i, 

4 = (4)^[l + 0.50226€a2]. 

(3.29) 


3.3 Near-Resonance Response Neglecting Longitu- 
dinal Inertia 

In this section, the near-resonance response of a harmonically excited non-linear beam 
is derived using the non-linear normal modes. The equation of motion of such a beam 
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can be written in the following non-dimensional form : 


d^w d^w 
dr^ dx^ 



= f{x) cos Or. 


(3.30) 


The response can be written as 

w{x, r) = 51 cos Or, (3.31) 

j=i 

where ^j(x) is the j-th non-linear normal mode, derived in the previous section. One 
should not assume that the non-linear normal modes are orthogonal to each other. 
When the forcing frequency is close to a linear natural frequency, say, wj,, then the 
participation of the corresponding non-linear normal mode is largest whereas the other 
modes participate only weakly. Mathematically this can be written in the following 
form : 

when Q = ojI, + €Qi, aj = 0{e) = ebj, (say), for j # n. (3.32) 


Substituting equations (3.31) and (3.32) into equation (3.30) and neglecting terms 
O(c^), one obtains 


CLn \ ^ '^n 


d^\ 

dx^ / 


cos — eal 


Jo ^ dx^ 


dx 


dtpn 


+e ^ bj 

ji^n 




A% 


dx^ 


dx^ 

cos flr = f cos f2r. 


cos^ Qr 


(3.33) 


Applying the harmonic balance technique to equation (3.33) and substituting equation 
(3.3), the following result is obtained: 

r 1 

(3.34) 


j¥^n 




= /■ 


Further, it can be noted from equation (3.13) that 

^ dx = 0(e), (j n). 


(3.35) 
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Hence, multiplying equation (3.34) by and integrating over x, while neglecting 
terms 0 (€^), one obtains 


(a;2 


(3.36) 


Equations (3.14) and (3.36) yield a cubic equation for the determination of a„. 

Since the non-linear frequency is close to the linear natural frequency and the n-th 
order linear mode of the beam is resonantly excited, one finds = 0(e). Now 

multiplying equation (3.34) by 'ipm, (tu ^ n), and integrating over x, one gets the 
following result using equations (3.35) and (3.32) and neglecting terms O(e^): 


fo f 


(m 7 ^ n). 


(3.37) 


It is observed from equation (3.37) that the same expression for Um (m # n) could 
have been obtained using only the linear theory. Substituting equations (3.36) and 
( 3 . 37 ) into equation (3.31), one finally obtains the near-resonance response. 

Now a special situation when some non-linear normal modes are orthogonal to each 
other is discussed. This happens, for example, for all the non-linear modes of a simply- 
supported beam (see Section 3 . 2 . 2 ) and the non-linear modes of even and odd orders 
for a clamped-clamped beam (refer Section 3.2.3). In such cases, when two orthogonal 
non-linear normal modes, say of order n and m are excited, the response can be written 
as 


with 


w{x, t) = Onlpn COS fir -l- am'fpm COS fir bjlpj cos fir, (3.38) 

j^riym 


[ ipni’m dx = 0 , [ Jpnfpj dx = 0 (e), / dx = 0 (e), when j 7 ^ m, n. 

^0 Jo Jo 

(3.39) 

Equation (3.38) is first substituted into equation (3.30). Then considering the fact that 
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for the cases under discussion, fg dx = 0 (for one obtains 

“ feOnC^ f dJ ^ 

4 Jo dx 


dx^ 




Jo^ dx’ 


^ j2 


dx 


dx^ 




d% 


j^m^n 


dx'^ 


= /, 


(3.40) 


by balancing of harmonics. Multiplying equation (3.40) by ipn and tpm separately and 
integrating over x, the following equations are obtained by neglecting terms O(e^): 

2 IT/, \ 2 


.^n ~ "I" 


^ /•(^) 0- 


Si't’lix 


Si Ax 

li^Ax 


(3.41) 


^m(^rn 


, [/o [ty H [j-o 


fo 


_ fJf(x)'l/Smdx 

foi’ldx ■ ^ - ’ 

Equations (3.41) and (3.42) can be solved simultaneously to obtain a„ and am- It may 
be noted that even the linear theory can be used to calculate a/s when j ^ m, n. 


3.3.1 Numerical Results and Discussion 

The non-linear response of a harmonically excited clamped-clamped beam was obtained 
in reference [146] using the linear modes. The response was assumed to be of the form 

3 

w{x,r) = '^A^(f>j{x) cos Q,t, 
i=i 

and the coefficients i4^’s were determined by solving three simultaneous non-linear 
algebraic equations which were obtained using Galerkin’s technique. 

Prom the analysis presented in Section 3.3, one can see the efficacy of using non- 
linear normal modes for calculating the response of a non-linear system. While the 
harmonic balance method presented in reference [146] requires solution of three simul- 
taneous non-linear algebraic equations, the present method reduces the computation 
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to merely solving one cubic equation. Furthermore, in the ordinary harmonic balance 
method, the number of equations increases proportionately with the number of linear 
modes considered. But in the method presented here, all the modal participations can 
be obtained from a single equation. 

To illustrate further that both the methods give same results, the near-resonance 
response of a clamped-clamped beam to a point harmonic excitation are calculated. 
The beam parameters are taken to be the same as those in reference [146]. To compare 
the results with those reported in the same reference, the coefficients of the first three 
linear normal modes (denoted by, say, Aj, with j = 1, 2, 3) are collected from equations 
(3.31), (3.36) and (3.37). With a forcing frequency close to the n-th linear natural 
frequency, the participation of different linear modes are as follows: 

An = an, ^;, = o„i + eAW, (m^n). 

For the dimensions of the beam treated in reference [146], the non-dimensional 
parameters are obtained as 7^ = 2.12 x 10~® and f{x) = 55624.77755 (a; — Xq), where 
xo is the point of application of the load. When the load is symmetrical about the mid- 
point of the beam and V’n is antisymmetric, the response amplitude is very small (since 
Gji = 0) and can be calculated from the linear theory. The situation is exactly similar 
when the load is antisymmetric and is symmetric. For the simply-supported and 
clamped-clamped beams, it was seen that V’n is symmetric or antisymmetric according 
to whether n is odd or even, respectively. Thus, from the above argument the non- 
linear analysis is required only if: (i) n is odd and the load is symmetric about the 
mid-point i.e., x = 1/2; (ii) n is even and the load is antisymmetric about a; = 1/2; 
(iii) the load is asymmetric. 

Figure 3.1(a) and 3.1(b) show the participation of the first and third linear modes 
when the first linear mode is resonantly excited by a load at a: = 1/2. The factor (7^), 
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Figure 3.1: (a) Participation of the first mode in the response, (b) Participation of the 
third mode in the response. A: linear theory, □ : Ref. [146] , o : present method; 
: in phase, : out of phase. 
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Figure 3.3: (a) Participation of the first mode in the response, (b) Participation of the 

third mode in the response. □ : Ref. [146], o : present method; : in phase, — 

- : out of phase. 
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used in the ordinate to compare the results of reference [146] arises because of the 
difference in the non-dimensionalization procedure. The results of the present analysis 
are in very good agreement with those of reference [146]. Both the results show 
marked deviation from the linear theory. It is also to be noted that the experimental 
data obtained in reference [146] support the theoretical results. Further, like all cases 
of cubic non-linearity, the intermediate value of the multiple solutions for the periodic 
response is unstable. 

When the first mode is resonantly excited by a point harmonic load applied at 
X = 1/4, both the first and second non-linear normal modes, which are orthogonal 
to each other, get excited. The response amplitudes in these modes are obtained by 
solving equations (3.41) and (3.42) and are shown in Figures 3.2(a) and 3.2(b). The 
results again show excellent agreement with those reported in reference [146]. 

When the third linear mode is excited by a point harmonic load at rr = 1/2, the 
first mode is also excited due to the modal coupling. The modal participation of the 
first and third linear modes are shown in Figures 3.3(a) and 3.3(b), respectively. The 
results also show reasonable agreement with those of reference [146]. 


3.4 Non-Linear Normal Modes Including Longitu- 
dinal Inertia 

In this section, the non-linear normal modes and the associated natural frequencies of 
an axially stationary beam are obtained after including the longitudinal inertia. The 
non-linear normal modes are then used to obtain the near-resonance response in the 
next section. The equations of motion of the beam, obtained by substituting c = 0 
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and To = 0 in equations (2.7) and (2.8), are 


1 d^u dw d'^w 

dr^ 2e dx^ ^ dx dx^ ’ 


(3.43) 


d'^w d^w 5 [ 1 fdu dw 


L/ uy \X/ V/ J. ^ V > V CL/ V ^L-ZCL/vJ 

dx^ dx 2e 5a; ^ dx ’ 


(3.44) 


where e = 7^/2. For obtaining the non-linear normal modes, the response w{x^t) 
assumed to be given by equation (3.2) is substituted into equation (3.43) to yield 


d^u 1 d'^u 2<^V’nd^^„ , o, . _ 

575 - “ IT'd?' + “ IT 

The steady-state solution to equation (3.45) is assumed in the form : 


(3.45) 


u{x, t) = ui{x) -f U 2 {x, r), 


(3.46) 


where ui(x) and U 2 {x,t) are the solutions with only the first and the second term, 
respectively, on the right hand of equation (3.45). Thus, 




(3.47) 


with g as the constant of integration. Equation (3.47), on integration and after using 
the boundary conditions that the longitudinal displacement is zero at both a: = 0 and 


z = 1, yields 


5 = € 


■a^f 


Next, the time-dependent part of u, U 2 {x,t) is easily obtained as 

. N 2^ CjTjix) . 


where 




(3.48) 


(3.49) 


(3.50) 
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and Fj’s and u/s are the linear mode shapes and natural frequencies for the longitudinal 
vibration, i.e., 


rj(x) = sin jnx and uj = jn, . j = 1, 2, 3, 


(3.51) 


Substitution of equations (3.2), (3.46)-(3.51) into equation (3.44) yields the follow- 
ing result after balancing the harmonic: 




dx^ 


= a 


+-a^ 

2 


dx 


4 I dx ' 


a 





2 


Jo ^ dx^ 


dx 


dx^ 


(3.52) 


For small non-linearities, the non-linear mode "tpn S'lid the associated frequency a;„ are 
expanded as in equations (3.4) and (3.5), respectively. Then substituting these into 
equation (3.52) and subsequently equating the coeflBcients of the like powers of e from 
both sides of the resulting equation, one obtains 


dVn 


0 f I \2± I ^ n 

€ • -(‘^n) (l>n + = 0, 


1 ! l\2s:(l) QS 2^-^ 1 2 


where 


A = 


1 

4 I dx , 


1 ^ ^3 
^J=l 


and 


Jo ^ dx ' 


dx 


d <f>n 

dx^ ’ 





(3.53) 

(3.54) 

(3.55) 


fl d^a d, <t>n p . 

dx (n CRN 

The solutions to the above equations are carried, out by the method described in 
Section 3.2. The final results are obtained as 


i’n — <j>n + ^ 




(3.57) 
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and 




(3.58) 


(3.59) 


where cnj^^’s are given by equation (3.10) and 

Now we recalculate the non-linear normal modes and the associated frequencies from 
equations (3.58) and (3.59) for simply-supported and clamped-clamped end conditions 
to demonstrate the effect of longitudinal inertia. 


3.4.1 Simply-Supported Beam 

Substitution of equation (3.22), i.e., the linear normal mode shape and frequency into 
equations (3.58) and (3.59) yields the following results 

1 + , ^n{x) = <i>n{x), (3.60) 

L o J 

which are identical with those given by equation (3.23). 


3.4.2 Clamped-Clamped Beam 

Using the linear mode shapes and frequencies from equation (3.24), the first three 
non-linear modes and the corresponding frequencies, are obtained as 

+ 0.006598€aV3, = (^i)^ [l + 0.235l€a2] , (3.61) 

-02 = 02j ^2 ^ (^ 2 )^ 0-41696(1^1 , (3.62) 

-03 = 03 — O.O5O62€a^0i, = (^^ 3 )^ [l + 0.502256a^j . (3.63) 

Comparing equations (3.61)-(3.63) with equations (3.27)-(3.29), one can conclude that 
the effect of the longitudinal inertia on the non-linear normal modes and frequencies 
of transverse vibration is negligible. 
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3.5 Near-Resonance Response Including Longitu- 
dinal Inertia 


If the beam is excited by a transverse harmonic force, the equations of motion are given 
by equation (3.43) and 


d^w d^w _ d 
dr^ dx^ dx 


1 .du.fdw. fdw^o 


+ f{x) cos Qt. (3.64) 

When + eOi, the response, as in Section 3.3, is written in the form: 


w{x, t) = anipn cos fir + e ^ bj'ipj(x) cos Qr. 

jjin 


(3.65) 


Equation (3.65) is first substituted into equation (3.43) and solved for u{x, r). The ex- 
pression for u(x, r) thus obtained is then used in equation (3.64). Thereafter, balancing 
the harmonics and omitting terms O(e^), the following equation is obtained: 


a„ -f2>n + 


. d 

dx^ ) ” da: 


^ ^ as 

4 da; y 2^^ u] 


£ (dil^nV 



2“" 




dx 


d^i>n 

dx^ 


+ €^bj 


j^n 




Substituting equation (3.52) into equation (3.66) one gets 


an^n(w^ - fi^) -H € ^ bj 

j^n 


-|- 


dVj 


dx^ 


dx^ 


= /• 


= /• 


Using equation (3.35) one finally obtains from equation (3.67) 

a fof(^)^ndx 

(^n-^^)fo^ndx' 

As noted in Section 3.3, 6j’s are obtained from linear analysis with j ^ n. 


(3.66) 


(3.67) 


(3.68) 
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3.5.1 Numerical Results and Discussion 


Here the results, reported iU' Section 3.3, are recalculated taking effects of the longitu- 
dinal inertia into consideration. The results are then compared with those of Section 
3.3. 

Figures 3.4(a) and 3.4(b) show the participations of the first and the third linear 
modes when the first linear mode is excited by a point load at a: = 1/2. The effect of 
longitudinal inertia is seen to be extremely small. 

Figures 3.5(a) and 3.5(b) show the participations of the first and the third linear 
modes when the third linear mode is excited by the same loading. Here also the effect 
of longitudinal inertia is seen to be negligibly small. 
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Figure 3.4: (a) Participation of the first mode in the response, (b) Participation of the 
third mode in the response, : in-phase response neglecting longitudinal inertia, - 

. out-of-phase response neglecting longitudinal inertia; o : including longitudinal 
inertia. 
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(a) 



/ 2 /®' 

Figure 3.5: (a) Participation of the first mode in the response, (b) Participation of the 

third mode in the response, : in-phase response neglecting longitudinal inertia, - 

— : out-of-phase response neglecting longitudinal inertia; o : including longitudinal 
inertia. 
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Chapter 4 


NON-LINEAR VIBRATION OF A 
TRAVELLING BEAM 

4.1 Introduction 

In this chapter, the concept of non-linear normal modes, explained for an axially sta- 
tionary beam in Chapter 3, is extended to the case of a travelling beam where even 
a standing linear mode is not possible. However, a complex normal mode can still be 
derived. Obviously, standing mode is not expected if the non-linearities are included in 
a travelling beam. Consequently, the existence of ‘non-linear complex normal modes’ 
are looked for. Two methods, namely, the modal analysis and wave-propagation ap- 
proach, are outlined in Section 4.2 to formulate these non-linear complex normal modes. 
Near-resonance steady-state response to an external harmonic excitation is presented 
in Section 4.3. Towards this end, the non-linear complex normal modes are used. As 
expected, the non-linear response has multiple solutions in some ranges of excitation 
frequency. A stability analysis, required to delineate the stable/unstable solutions, is 
also carried out in Section 4.3. 
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For a near-resonant excitation, the contribution of one linear mode is order of mag- 
nitude higher than the others and the modal analysis converges quickly. However, for a 
hard non-resonant excitation, a number of linear modes participate with significant con- 
tributions and the modal approach becomes unwieldy. Such an excitation is considered 
in Section 4.4. The problem of modal coupling is alleviated by taking the linear re- 
sponse in a closed form as obtained in Chapter 2 with the help of the wave-propagation 
theory. Thereafter, the non-linear response is obtained using a perturbation approach. 


4.2 Complex Normal Modes 


In this section, the non-linear complex normal modes and the associated non-linear 
natural frequencies of a simply-supported travelling beam are derived. The equation 
of motion of the non-linear free vibration is given by equation (2.16) and the boundary 
conditions by equations (2.10) and (2.11). 

To analyse the non-linear free vibration, one can write 

W{x, t) = ^^n(a;)e“’*^ -H c.c. (4.1) 


where Un is the non-linear natural frequency, and = 


ioJni’n 


1 , with rpn as the non- 


\ ] 

linear mode shape, and c.c. denotes the complex conjugate of the preceeding expression. 
Both u)n and xpn are dependent on the complex amplitude a. 

Substituting equation (4.1) into the equation of motion (i.e., equation (2.16)) and 
balancing the first harmonic, (or equating the coefficients of and separately, 

to zero), one obtains the following complex equation and its conjugate: 

ia;„|A^n-b^B^„-€Ni = 0, (4.2) 
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where 


Nx 


j a^a ( f^Vdx 
’ " ' dx drc ^ ^ dx 2 Jo V dx J 


8 


dx^ 


■•)■ 


(4.3) 


For weak non-linearity, it can be assumed that the non-linear modes are very close to 
the linear modes and one can write 


+ • • • , 


(4.4) 


and 

uJn = i^l + + • • • ■ (4-5) 

Substituting equations (4.4) and (4.5) in equation (4.2) and equating the coefficients 
of like powers of e from both sides, one gets 


£» : = 0 , 

+ |ba 1"> = -< ji3j”>A4„ + Ni', 


(4.6) 

(4.7) 


where 


Ni^ = 


’"/o’ 


8 I dx^ Jo dx dx dx^ Jo \ 




(4.8) 


Equation (4.6) is trivially satisfied. To solve equation (4.7), one can expand in 
terms of the linear modes as 


^(n) _ ^ + X) (4.9) 

m^n n=l 

since the solvability condition of equation (4.7) demands that A^"^ can not have any 
contribution from 

The right hand side of equation (4.7) is first expanded in terms of the linear complex 
normal modes (see Section 2.3.1). Then substituting equation (4.9) into equation (4.7) 
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and equating the coefficients of and ^j, {j = 1, 2, 3, • • •) from both sides of the 
same, one finally gets 


' ai/o'CA^ndx’ 
g’^ = — 


(4.10) 


(4.11) 


-, m= 1,2,3, 


(4.12) 


= 1 , „ = 1, 2, 3, . . . (4.12) 

Expanding and Ni' using equations (2.19) and (4.8), respectively, equations (4.10)- 

(4.12) take the following forms: 


M A 4>n) 

aa 4 - < ;> <te 


(4.13) 


(4.14) 


where 


hm. — 


_9'm _ 1 K^rn , <f>n) 


aa 4wJi + a;f,jj-^i$^A$TOdx 


, m = l,2,3. 


(4.15) 




Finally, equations (4.13)-(4.15) can be used to write the non-linear normal modes and 
the associated frequencies as 


+ eaa j 53 9m^m + X) ^m^rn j + 0{e^) 

\m^n m=l / 


(4.16) 




(4.17) 
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It can be easily shown from equation (4.13) that for a stationary beam 



^hSo4>l^x 


yielding 



a result obtained explicitly in Chapter 3 (see equation (3.20)) with the amplitude a as 
real. 

As in the case of an axially stationary beam, the non-linear complex normal modes 
of a travelling beam can be derived also by the wave-propagation analysis as presented 
in the next section. 


4.2.1 Wave-Propagation Approach 

During vibration in any of the non-linear normal modes, the phases of the waves 
corresponding to all the linear normal modes get closed. Assuming the non-linear 
complex modal response as : 

W(x, t) = (4.18) 

the following equation is satisfied during the modal vibration : 

= 0. (4.19) 

where Ni is given by equation (4.3). 

FVom the linear analysis presented in Chapter 2, it is observed that to close the 
phases of the waves associated with and in time 27r/f2, the required forces are 

{/} = i{Q, - 
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and 


respectively. As the force due to the non-linearity closes the phases of the waves of all 
the linear modes simultaneously, one can write, 

oo ^ oo 7 

eNi = i(u‘^ - a;)-^A^rn - '£ (4.20) 

m=l Tn=l ^ 

As in the previous section, for weak non-linearity, one assumes 

u; = col + 0(e) (4.21) 


and 


m=l m=l 


(4.22) 


where a„ = a, am = 0(e) for m n and bm = 0(e) for all m. Combining equations 
(4.20)-(4.22) and the orthogonality relations, one gets the following results, valid up to 

o(e) : _ 

7„'#jNldi 


and 


a. 


bm — 


u 


2e 






/o'CN'.dc 


- “m) /o da: 

2e /„■ «^N! dx_ 


m^n 


i (w!. + a,<.);„'4^A«„da:’ 


m = 1,2,3, 


where is the vector Ni with changed to 0„. It can be verified, by expanding 
$j’s and that the above results are identical to those obtained in Section 4.2. 


4.3 Near-Resonance Response 

In this section, the non-linear response of a travelling beam, excited harmonically with 
a frequency very close to one of the linear natural frequencies, is detailed. 
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The equation of motion (i.e., equation (2.17) is reproduced below : 


A— + BW-eN = i, (4.23) 

where 

f = {f{x) cos fir , 0}^ . (4.24) 

Equation (4.24) can be written as 

f = + c.c. (4.25) 

with fi = {f{x)/2 , 0}^. 

The principal harmonic response of the beam can be assumed as 

W(x, r) = i (A(rc)e'“^ + A(x)e-'^") . (4.26) 

Now, in the absence of any internal resonance, if the excitation frequency is very close 
to a linear natural frequency, then the response contains primarily the corresponding 
non-linear mode. The contributions of the other non-linear normal modes are very 
weak. This statement can be expressed mathematically as follows: 

When fl = -t- efli, 

A = Oji^n + ^ + 0(e^), (4.27) 

m^n m=l 

with a'^ = 0(e) = eOm, (say) and = 0(e) = ebm, (say). 

Equations (4.26) and (4.27) are then substituted into equation (4.23) and the har- 
monic balance method (i.e., equate the coefficients of and separately, from 
both sides of the equation) is used. Neglecting terms O(e^) the following relation and 
its complex conjugate are obtained: 


-f- 

rnsin ^ ^ I J 
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(4.28) 


+ E 

m-l \ ^ ^ } 


= fi, 


where Ni is defined as 


(#r sf- S'/; (S)’ *•! 4’- 


V 

Substitution of equation (4.2) in equation (4.28) yields 


i (f2— CUn) A^n+C 




[m^n 


2 


)+e( 


. * -3- . 




(: 


= fl. 
.29) 


Premultiplying equation (4.29) by and using equation (4.16) and the orthogonality 
relations (see equations (2.27) and (2.28)), one gets, after neglecting terms O(e^), 


CLn — 


2/^^;;;fidx 


(4.30) 


i(^^-^n)/o<A^„da;' 

Since ujn is related to a„ through equation (4.17), equation (4.30) is a complex cubic 
equation in a„. To solve equation (4.30) one can write 


ttfi — u^e “. 


Now, expanding and fi with the help of equations (4.13)-(4.17) and (4.25) and 
again considering terms up to 0 (e), the following results are obtained: 

Q + ealM 


On (cos 6n + i Sin dn) = 


{Q-ul-ealS)tn' 


(4.31) 


where 



and 


tn= f ^n-^^ndx. 

JO 

It is to be observed that Q and M are complex quantities whereas S and are real. 
So, equating the real and imaginary parts of equation (4.31), one gets 

+ ealM^ 


and 


CLfi COS Oji 


a„ sin On = 


- ojJj - ealS)tn 
+ ealM^ 


(4.32) 


(4.33) 


(JI - w', - €alS)tn ’ 

where the superscripts R and I denote the real and imaginary parts, respectively. 
Squaring equations (4.32) and (4.33) and adding, one gets, 

(Q^ + ea^M^y + (Q^ + ea^M^y 

Again equation (4.34) is a cubic equation in o^. If the externally applied force f(x) is 
a point load expressed by Fo5{x — a;o), equation (4.34) can be written as 


a? 


(4.34) 








Kril-rif 


S2 




2u{MfQf + M[Ql) 
(4.35) 


where 


n = J- = (0 >1 = -M., Ml = M/Fo, Qi = Q/F^. 

The phase angle On can be obtained from equations (4.32) and (4.33) as 

(Q^ + ealM^) 


tan On = 


(Q* + ea^M^) ‘ 


(4.36) 


Now, to solve for a'^ (appearing in equation (4.27)) from equation (4.29), one premul- 


tiplies equation (4.29) by and uses equations (2.27) and (2.28) to get 

m^n. 


CCLfn — 


2/o<fidi 




(4.37) 
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In a similar fashion 6^ can be obtained as 


brr, 




Si dx 


m = 1, 2 ■ 


(4.38) 


Comparing equations (4.37) and (4.38) with equations (2.57) and (2.58), respectively, 
one can see that the same expressions for 5^ and a'^ could be obtained using only the 
linear theory. 

Collecting On, 6n, o!m from equations (4.34)-(4.38) and putting them into 

equations (4.26) and (4.27), one gets the non-linear response of a travelling beam to 
a near-resonance harmonic excitation. Obviously, the linear response, presented in 
Section 2.4, follows by subtituting e = 0. 

Since a continuous phase difference exists in a travelling beam, different sections 
reach the maximum amplitude at different instants. Hence, a better representation of 
the response is to show the envelope which bounds the response at all times. When the 
n-th complex normal mode is excited, the shape of the beam can very well be expressed 
in the form 

w{x, r) = -I- c.c. -f- (9(e) 

= a„0;cos(Or-h0n + e) + O(e), (4-39) 


where 


and 

<f>i = 


4.3.1 Stability Analysis of Periodic Solutions 

As expected for a non-linear system, the possibility of a multi-valued response to a 
harmonic excitation is observed in Section 4.3. This is manifested as a cubic equation 
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where Li = {li(x,r) ,0}^, with 


2i(nr+&„) 

dx2 


+ -a^ 


l_o \dVn 


4“" 1 \Jo ^ dx ^ dx^ 


+ 2 f /' dx) 

\^o dx dx J dx2 J 


+ 




dVn ' 

dx2 


Applying the orthogonality conditions and equation (2.20) in equation (4.43), one gets 
the following equation and its complex conjugate: 


d^n . If. , Jo dx 


/o^$„Lidx - 


d'T /o dx /o dx 


^n- 


(4.44) 


Assuming r' = Qr + 6n and equation (4.44) can be written as 


O 


dr' 


-iV^n 




Si «„L'i di - 


C I ^ 

7„»<A*„da: ” ;„*CA4„d2: 


^Tl) 


(4.45) 


where L'l = {Zi(x, r') ,0}^. 

Since = o;^ + 0(e), onq can assume 




(4.46) 


and 


77 = 1 + e 77 i 4 . (4.47) 

Substituting equations (4.46) and (4.47) into equation (4.45) and equating the like 
powers of e from both sides, the following results are obtained: 


dr' 


(4.48) 


dr' /o'$„A$„dx 


^nL'i dx ^ 0 ) 


(4.49) 
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or 

J7i = ri-r2. (4.53) 

Putting J? = ^ and equations (4.52) and (4.53) in equation (4.47), one gets the bound- 
ing frequencies as 

fi' = wi-£(r,+r2K (4.64) 

and 

n" = u>i ^ e(ri - r2)wj,. (4.55) 

When Pi and r 2 are substituted from equation (4.51), the following relation is finally 
yielded: 
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where is given by equation (4.13). This implies that equation (4.55) is nothing 
but the backbone curve. The primary unstable region is bounded by the two curves, 
expressed by equations (4.54) and (4.55), in the Q-a„ plane. 

4.3.2 Numerical Results and Discussion 

In this section, some numerical results obtained from the theoretical analysis presented 
in Sections 4.2 and 4.3 are discussed. The initial tension, To is taken to be unity. 

Figure 4.1 shows the amplitude dependence of the first non-linear frequency, ex- 
pressed as a factor of the first linear natural frequency of a stationary beam (w®). 
The non-linear frequency is obtained using equations (4.13) and (4.17) with = aa. 
The figure shows a hardening type non-linear characteristic. Consequently, the non- 
linearity, as expected, has a stabilizing effect on the divergence instability. The hard- 
ening effect becomes more pronounced with increasing travelling speed of the beam. 

Using the non-linear normal modes, discussed in Section 4.2, the steady-state re- 
sponse of a travelling beam, excited by a point harmonic load applied at xq = 1/3, 
and having a frequency close to the first linear natural frequency (i.e. n = 1), is calcu- 
lated taking the effects of the first two linear modes. The roots of A (i.e. from 

equation (4.35), for three values of the travelling speed d = c/(ccr)i = 0, 0.3 and 0.5, 
with Foy/e = 13.22, are shown in Figure 4.2. The frequency of excitation, fi, has been 
expressed as a factor of the first natural frequency of a stationary beam ( 0 ;°). Since 
the mode shape of a travelling beam is not stationary in space, the plot of only di 
becomes meaningless. In order to provide useful information about the response, the 
positive half of the shape of the response envelope, i.e., ^^(a:), shown in Figure 4.3, is 
also necessary. It is to be noted that the envelope bounds the beam displacement 
within an accuracy of o(l), and different sections of the beam touch the envelope at 
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different instants of time. Hence, Figures 4.2 and 4.3 are to be used together to obtain 
the response envelope for different travelling speeds and excitation frequencies. 

The results of the stability analysis of the steady-state harmonic solution, presented 
in Section 4.3.1, is shown in Figure 4.4. The frequency response curves for a beam 
having c = 0.5(Ccr)i and Fo\/c = 12 and 53 are plotted. The intermediate branches of 
both the plots, shown by the broken lines, refer to unstable solutions. As explained in 
Section 4.3.1, the region between the curves, represented by equations (4.54) and (4.55), 
refers to unstable solutions. These boundary curves are shown by chain lines. The 
upper boundary of the unstable region, shown in Figure 4.4, represents the backbone 
curve. The other curve is the locus of the points of vertical tangencies of the frequency 
response curves. 


4.4 Response to Non-Resonant Hard Excitation 

In this section, the non-linear stady-state periodic response of the beam subjected to a 
non-resonant hard excitation is analysed. As the frequency of excitation is away from 
any of the linear natural frequencies or combinations thereof, various linear modes are 
excited. Consequently, the modal analysis is not convenient for obtaining the response. 
In such a situation, the transfer function of the beam, derived in Chapter 2, can be 
used with great advantage. Since the linear response in the present method, contrary 
to the modal analysis, is obtained in a closed form, the analysis becomes simplified. 
Also the error encountered in the modal analysis due to neglecting the higher order 
modes is avoided. 

Let the excitation be of the form f{x, r) = FoS{x — rro) cos fir, with Fq sufficiently 
large so as to render the effects of the non-linear term in equation (2.15) to become sig- 
nificant, even if the excitation frequency Q, is away from any of the natural frequencies 
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Figure 4.4: Primary unstable region in the frequency-response diagram of a resonantly 
excited non-linear travelling beam. To = 1.0, c = 0.5(ccr)i. I : F^y/e = 12, II : 
Fq^/c — 53. : stable solutions, : unstable solutions. 
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or their combinations. The response of w{x,t) is sought in the following series form : 


w(x,t) = Wq{x,t) + ewi{x,T) + • • • . 


(4.56) 


It may be pointed out that the above series is valid only for Q away from any linear 
natural frequency, When the series diverges. 

Substituting equation (4.56) into equation (2.15) and equating the coefficients of 
the like powers of e from both the sides, the following results are obtained : 


c»: ^+2c 


A . 


dr^ 


d^wp 

dxdr 


d^wp d^wp 


+ (c -ro)^ + 


dx^ dx^ 


= f{x) cos fir 


d^wi d^wi 


+ 2c^ + (c^-To):^ + 


^,dwp 

dx 




)^da; 


d^wp 

dx^ 


(4.57) 


(4.58) 


dxdr ' dx"^ ' dx^ 

Following the procedure outlined in Section 2.4.2, equation (4.57) is solved to yield 

wp{x,t) = — Gix, C, iQ)f{OdC + —^ G{x, C, -in)/(C)dC (4.59) 


or 


Wp{x, t) 


H{x, i^) mr H{x, -iiiT 


(4.60) 


where 

H{x, iCl)== fG{x, C, if^)/(C)<iC- 
Jo 

Solution of equation (4.58) can be obtained by writing the right hand side of the same 
equation as follows : 


where 



= [Mi(rr)e*^^ + M2e3*'^^] 


+ C.C. 


(4.61) 


A/,(x) = - 


f 


fdH{x, iQ) 
V dx 



d^H{x, —ifl) 
dx^ 



dH(x, iQ) dH(x, -ifi) 
dx dx 


d^Hjx, iQ) 
^ dx^ 
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and 


M2{x) = 


8 


d‘^H{x, in) 


dx"^ 


With the help of equation (4.61), Wi{x,t) is obtained as : 

= G{x, C, if2)Mi(C)dC + e3‘^^ G{x, C, 3ifi)M2(C)dC + c.c.. (4.62) 
Jo Jo 

Substituting equations (4.60) and (4.62) in equation (4.56), the response of a non-linear 
beam, up to the term o(e), is obtained. Further simplification can be obtained if the 
force is assumed to be a point harmonic one i.e., f{x) = FoS(x — Xq). In that case 


H{x,in) = FqG{x, Xq, in), 


and the values of M\{x) and M^ix) are evaluated by numerically computing a few 
definite integrals. 


4.4.1 Numerical Results and Discussion 

The non-linear response of a beam subjected to a point harmonic excitation is presented 
in this section. The axial speed, initial tension, amplitude of excitation and the small 
parameter e are taken as follows : 

c'=0.5. To = 1.0, ^0 = 1000, e = 10-2. 

The Fourier transform of the response at the point a: = 0.75 to a non-resonant hard 
excitation applied at xq = 0.3 is shown in Figure 4.5. The response has both harmonic 
and superharmonic components. Figure 4.6 shows the variations of harmonic compo- 
nents {w^i{x,t), say) of the response amplitudes, measured at two different locations 
(x = 0.5 and x = 0.75) of the beam, to a point excitation having away from both uj[ 
^2 (=38.85552). It is seen that the effect of the non-linearity is more 
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3.0 



Q 


Figure 4.6: Variation of the harmonic components of responses to point harmonic 

excitation. Fq = 1000, (/ = 0.5, e = 0.01, xq = 0.3; : non-linear response 

measured at x = 0.75; — ; linear response measured at x = 0.75; ; non-linear 

response measured at x = 0.5; ■ — ; linear response measured at x = 0.5. 
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pronounced at x = 0.75, since the participation of the second mode is stronger for that 
location. 

As discussed earlier, the non-linear response can be obtained by either the wave- 
propagation method or the Galerkin’s technique. Figure 4.7 shows the amplitude of the 
harmonic component of the response calculated by these two methods. The discrepancy 
between the results obtained by these two methods is seen to be more than that for a 
linear beam. This is perhaps due to the coupling between various linear modes in the 
presence of non-linearity. 
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Chapter 5 


PARAMETRICALLY EXCITED 
TRAVELLING BEAM WITH 
AND WITHOUT EXTERNAL 
FORCING 

5.1 Introduction 

In the last chapter, the effects of the non-linear term on the forced response of a 
travelling beam were analysed. In this chapter, another form of excitation, namely, the 
parametric excitation, is considered. 

The axial movement of a travelling beam is normally sustained by rotating pul- 
leys which are either rigidly or flexibly mounted upon the foundation. Normally, a 
finite pulley-support compliance in the axial direction provides a stabilizing effect to 
the buckling or divergence instability that occurs in the super-critical speed regime. 
However, the presence of any rotating unbalance in the flexibly mounted pulleys may 
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cause parametric excitation, which in turn causes instability in an otherwise stable sys- 
tem. As discussed in Chapter 1, there may be other causes of parametric excitations 
in various travelling slender members. 

Non-linearities have long been recognized to play a significant role in the para- 
metrically excited Duffing oscillator [145]. Whereas the linear analysis predicts an 
exponential growth of the amplitude of such an oscillator for a range of excitation am- 
plitudes and frequencies, the non-linear terms limit the response to a finite amplitude in 
the steady-state. These solutions are referred to as ‘limit cycles’ in the literature. The 
existence of limit cycles in a travelling string has already been mentioned in Chapter 1. 

The response of a Duffing oscillator to simultaneous external and parametric ex- 
citations [147, 148] is normally aperiodic. But if the frequency of the two forms of 
excitation are synchronized, a steady-state harmonic solution is obtained. Superhar- 
monics and subharmonics may also result if some special relationships exist between 
the frequencies of the two forms of excitation. 

After deriving the relevant equations of motion in Section 5.2, the primary resonance 
of a travelling beam under only parametric excitation including the non-linear effects is 
presented in Section 5.3. To this end, the concept of non-linear complex normal modes, 
explained in the previous chapter, has been used. The limit-cycle amplitudes and their 
stabilities have been given special attention. The excitation of other linear modes 
besides the primarily excited one has been clearly brought out. Numerical solutions 
are presented to show the effects of different system parameters (like band-speed, pulley 
compliance etc.) on the limit-cycle amplitudes. In Section 5.4, the non-linear complex 
normal modes are again used to obtain the steady-state harmonic response of the 
beam when both forms of excitation, namely parametric and external harmonic, are 
simultaneously active. 
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5.2 Equation of Motion 


Consider a travelling slender beam, driven like a belt-pulley system as shown in Figure 
5.1. The finite pulley compliance is modelled by a linear spring having stiffness A;*. 
Under the usual assumption u* = 0(w*}^, the equations of motion for undamped 
planar vibration including the non-linear terms are: 


pA 


^ ,d^u* ,2d^u* 

dt^ ^ d^dt ^ de 


EA 


d^u* 

W 


(EA-T*) 


dw* d^w* 

ae de 


(5.1) 


and 


pA 


d'^vj* 

~W 


“t* 2c 


,d^w* 

d^dt 


+ c* 


'.d^w* 


de 


-T 


d^w* 


de 


EL 


d^w* 


= {EA - T*) 




du* dw* 1 / dw* 


(5.2) 


where T* is the equilibrium tension in the beam and the other symbols used in equations 
(5.1) and (5.2) are explained in Chapter 2. The equilibrium tension T* and the initial 
tension Tq are related as [37] 


T* =T^ + r}pAc*\ (5.3) 


where 

Neglecting the curvature of the beam beyond the frictionless guides, the boundary 
conditions can be written as follows: 

w*{0,t) = w*{l,t)=:0, (5.5) 

d^w*{0,t) _ d^w*{l,t) _ p 

~~de de~ 

u*{0,t)=0 
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(5.6) 

(5.7) 




Figure 5.1: Schematic diagram of a travelling beam having parametric excitation. 
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and 


du* 1 



EA^ 




cos 2Q.pt, 


(5.8) 


where m*e\ is the unbalance of the pulley having a rotational speed 20*. Equation 
(5.8) is obtained by balancing the horizontal forces at the elastically mounted pulley. 

Attention may be drawn to the fact that in the absence of any frictionless guide, 
the dynamics of the two spans of the beam/wheel system are coupled [43, 44, 46, 47] 
owing to the finite curvature at the beam-wheel interface. The coupling becomes less 
significant with increasing wheel radius and/or static tension and also with decreasing 
flexural rigidity of the beam. It is well known that the static deformation becomes 
significant at an axial speed close to the critical value. However, this deformation is 
reduced by incorporating the frictionless guides. Consequently, the simply-supported 
boundary conditions for the transverse vibration can now be used without much error. 

The non-dimensional parameters used subsequently are given below: 




u = -, 7 = I = |, r = (^)‘/V/i, c = K = KII(2EA), . 




,e: 


Q;t = QpT, 2m* = m/ipAl), ei = el/{^% To ^n/iEA'f), R = R*/L (5.9) 


It may be pointed out that the relative magnitudes of non-dimensionalized parameters 
u and w are now related as u = 0{'y^w^). 

Using equation (5.3) and the non-dimensional parameters, the equations of motion 
and the boundary conditions are rewritten as 

d‘^u 2 ^^ 

^dxdr ^ dx^ 

dxjy 
(5.11) 

w(0, r) = w{l, r) = 0, (5-12) 


d'^w d^w 
+2c- 


dr^ dxdr 


+{kc^-To) 




1 du dw 'y'^ 
7 ^ dx dx 2 


1 d^u 
72 dx^ 


7^(1 - 'fn - 'f^n) 


dvj d‘w 
dx dx^ ’ 


(5.10) 
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t) _ 5 ^w(l,r) 
dx^ dx^ 


(5.13) 


^(0>’’) = 0. (5.14) 

k 

+ y“(l)7-) = 'r^mQlei cos 2f2pr, (5.15) 

where k = 1 - rj. It is easy to see from equation ( 5 . 4 ) that k = 1 implies a rigidly 

mounted pulley i.e., k* -> oo. Moreover, from equations ( 5 . 4 ) and ( 5 . 9 ) one notes that 
kp and k are related as 


du 7 ^ / dw \ ^ 
dx'^ 2 V^a:/ 


^p/ (^ + kp) — k and 1/ (1 + A^) = 1 — k. (5.16) 

Since 7 < 1 implies that the speed of the longitudinal waves is much higher than that 
of the transverse waves, one can neglect the longitudinal inertia term in equation ( 5 . 10 ) 
to get (neglecting terms of order higher than 

1 d^u ^ « dw d^w 

^ dx'd^^ ( 5 - 17 ) 

which when integrated twice results in 

= (^) *i+*Si(r)+« 2 {T). (5.18) 

gration constants gi(r) and q 2 {T) are obtained, using boundary conditions 
(5.14) and (5.15) in equation (5.18), as 


? 2 (r) = 0 , 


(5.19) 


and 




where eo = mQlei. 


(Tny ^ cos 2flpr, (5,20) 


gets 


Combining equations (5.11), (5.16), (5.18)-(5.20) and neglecting 


terms 0 ( 7 ^), one 


d^vj d^w 

^ ~ ~ 2 eeo(l - k) cos 2QpT)~ + ^ 
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where c(= 7 ^/ 2 ) will be treated as the small parameter, since c <C 1. 

The following points are at once observed from equation (5.21): 

a) As expected, the spring flexibility makes the system softer. 

b) For a rigidly mounted pulley (i.e., k = 1), the unbalance has no effect. This is 
because the unbalance force gets balanced by the reactive force at the pulley-support. 

A small viscous damping force efj.{dw/dT) is introduced into equation (5.21) which 
can then be recast in the standard state-space form as 


dW 

+ BW = eNi -H cNj, 

OT 


(5.22) 


where 


and 



O 


1 

1 


dw 

dr 

A = 


, B = 


, W = < 

f 


1 

o 

1 


1 

1 

o 


1 ^ J 


( d^w dw 

Ni = |2eo(l - k) cos 2npr^ - 




(5.23) 


(5.24) 


(/o (S) 

with / as the identity operator, K = {kc^ - ^o)|? + ^ 

In the presence of a harmonic excitation, expressed in the non-dimensional form as 

f{x) cos (fi/T + 6f), the equation of motion (5.22) becomes 

A^ BW = eNi + eN2 + f, 
dr 


f = {f{x) cos (fl/T + 9f), 0}^ . 


where 


(5.26) 



5.3 Limit-Cycle Amplitude without External Forc- 


ing 

Assuming that the amplitude of the solution to equation (5.22) changes slowly with 
time, one can write the response as 

1 


W(a:,r) = -Ape*”'-^ + c.c., 


(5.27) 


with An 




where Ap is a complex shape function. The abbreviation c.c. 


denotes the complex conjugate of the preceeding term. 

Depending upon the relation between fip and the linear natural frequencies (w^) 
of the travelling beam, different kinds of parametric excitation may appear. In what 
follows, only the primary excitation is analysed. 

To study the primary resonance, the following condition is first set 


Qp = Jn + ea„. 


(5.28) 


where c7-„ is called the detuning' parameter. In such a situation, the n-th mode is 
primarily excited and one can write 


Ap = a„(r)'If„ + € 


S) am(T)^,„ + ^ 6m (r)^' 

m^n rn=l 


m 


(5.29) 


where denotes the n-th non-linear complex modal vector with corresponding non- 
linear frequency a;„. Following the steps outlined in Section 4.2 one can show that 
and are given by equations (4.13)-(4.17) with a factor k in A(<^m, 4>n). Further, 
satisfies the following equation: 




(5.30) 
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where 




2t/i 


ka^a 

Y~ 


\Jo \dx J j dx^ yJo dx dx 


dx 


dlpn 

dx^ 


,0 


(5.31) 


For a slowly-varying amplitude 

daj(T) 


dr 


and 


d&j(-r) 

dr 


= 0(e) = eaj (say), for j = 1,2, 


0(e) = eSj (say), for i = 1,2, 


(5.32) 


Now, equations (5.27), (5.28) and (5.32) are substituted in equation (5.22). Thereafter, 
equating the coefficients of from both sides and neglecting terms O(e^), one gets 




m^n ^ 


where 


»n=l ^ 

aUl - 


eNi, 


(5.33) 


It is to be noted that while deriving equation (5.33), the fact that (f2 — a;„) is 0(e) (as 
can be easily seen from equations (5.29) and (4.17)), has been made use of. 

Using the orthogonality relations between the linear modes ($n) iJi equation (5.33), 

one gets 


. X 2e/o^$rNi'dx 

eOn + ian[}lp — (^n) — 7 > 


(5.34) 


where i„ is defined in Section 4.3. For notational.simplicity, o„ is replaced by a, and 
the complex amplitude a is expressed as 

(5^35) 


a = oe‘®, 


with a as real. 
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Expanding Wn as in equation (4.17) and using equations (4.13), (5.28), (5.32), and 
(5.35) in equation (5.34) and thereafter separating the real and imaginary parts, one 
gets 

— = — eenO \G^ cos 26 — sin 2d] — efj, a, (5.36) 

dr '• 

and 


a— = e——dn — ea^n + ’ 

dr 4 i„ '• 


(5.37) 


where 


dji “ k 


,Jq V dx ^ 


dx 


f d(f>n 

Jo \ dx 




lo dx dx J 


G = G^ + iG^ = ^^ dx, 

H= f Mn da;. 

JO 


and 


j _ eo^n(^ ~ 
in 

./ _ 

^ ~~Tr- 


(5.38) 


The steady-state solutions of equation (5.34) correspond to the fixed points of equations 
(5.36) and (5.37) and are obtained by putting 


da d9 

— = 0 and — = 0 (if a 7 ^ 0 ). 


(5.39) 


A ready inspection reveals that a = 0 is a fixed point, which is called the ‘trivial limit 
cycle’ in the literature [145]. The amplitude (a) and the phase angle (^) of non-trivial 
limit cycle(s) are obtained from equations (5.36)-(5.39). Thus, one finally obtains 


sin {29 — S) = ~ 




+ {O'f ’ 


(5,40) 
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cos {26 - 5) = — 


with 


+ (G')’ 


. = arctan(|l' 




4 


(6,41) 


Squaring and adding equations (5.40) and (5.41), one gets 

4t„ 


{ayftf 


f^ndn L 


e(T„ ± ^(€e'o)2 [(G^)2 + [G^f] - (e/i')2 


(5.42) 


Remembering that should be positive, it is easy to note the following three possi- 
bilities from equation (5.42) ; 

(a) No non-trivial limit cycle exists : 

(i) if 


eo < 




#¥+CT’ 

(ii) even if condition (5.43) is violated but < 0 and 


(5.43) 




(5.44) 


(b) Only one limit cycle exists, if condition (5.43) is violated and 

kJ < ^/et{G^Y + iG^r-t^'^. (5-45) 

(c) Two limit cycles exist, if both conditions (5.43) and (5.45) are violated for cr„ > 0. 
It must be emphasized that even if a non-trivial limit cycle exists, it may not be 

stable. The analysis detailed below is then required to confirm the stability of the limit 


cycle. 


5.3.1 Stability Analysis of Limit Cycles 

The stability of the limit cycle(s), obtained in Section 5.3, is determined by linearizing 
equation (6.34) about the fixed points of a and d. Since for the trivial Umit cycle 
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(a = 0) no unique phase angle (9) exists, the linear part of equation (5.34) is anlysed 
rather than using equations (5.36) and (5.37). However, for non-trivial limit cycles, 
linearization of equations (5.36) and (5.37) are carried out around a and 0, given by 
equations (5.40) and (5.41). In what follows, first the stability of the trivial limit cycle 
and then that of the non-trivial one are discussed, 

(1) TRIVIAL LIMIT CYCLE. 

Expand a as 

a = ax + iciy. (5.46) 

Now, equation (5.46) is substituted in equation (5.34) and is expanded as equation 
(4.16). Thereafter, neglecting the non-linear terms and equating the real and imaginary 
parts separately from both sides, one gets 

d [ 0:6 1 \ (C-D) (E - F) 1 f 0,6 1 . . 

[-{B + J’) -(C + D) J ( 0, J 

where C = e'^G^, D = , E = F = e'^G^. 

The eigenvalues of the coefficient matrix are then given by 

S 1.2 = -D± + [(C2 - L»2) - (F2 - F2)] , , 

which implies that the trivial limit cycle is unstable if 

K((GV + (GV)-m'‘]. (5.48) 

Comparing condition (5.48) with conditions (5.43)-(5.45), it is inferred that the trivial 
limit cycle is stable, if this is the only limit cycle or if two more non-trivial limit 
cycles exist. However, the trivial limit cycle is unstable, if there exists only one more 
non-trivial limit cycle, 

(2) NON-TRIVIAL LIMIT CYCLE. 

To ascertain the stability of the non-trivial limit cycles, one perturbs the amplitude 
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(a) and phase (0) in equations (5.36) and (5.37) in the form 


Up = a + r)p 


(5.49) 


and 

ep = e + Cp. (5.50) 

Substituting equations (5.49) and (5.50) into equations (5.36) and (5.37) and using 
equations (5.40)-(5.42), the following results are obtained (after linearization) : 


dr 


Vp 

i Cp 


= e 


0 2eoa^(G«)2 + (G^^ cos {29 - 5) 


1 ul.d„ 

I 2 tn 


a 


-2\l 


i 

Vp 

1 


Up J 


(5.51) 


(5.52) 


The eigenvalues of the coefficient matrix of equation (5.51) are 

= -ll ± + 4^5" KG*)'' + cos (2« - i) . 

FVom equation (5.52) one can show that the limit cycle is stable if cos (20 - 0) < 0. 
From equations (5.40) and (5.41) and conditions (5.43)-(5.45) the following conclusions 

are apparent: 

(i) When two non-trivial limit cycles exist, the value of cos (20 - S) is negative for 
the limit cycle with the larger amplitude and positive for the one with the smaller 
amplitude. Hence the limit cycle with the larger amplitude is stable and the other one 

unstable. 

(ii) When a single non-trivial limit cycle exists, it.is always stable becanse the valne ot 
cos (20 - 5) associated with it is negative. 

5.3.2 Numerical Results and Discussion 

Numerical results are presented to show the effects ot various system parameters on the 
steady-state solution, i.e., the amplitude of the limit cycle. AU the results have been 
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obtained with To = 1. In what follows, the effects of only those parameters which can 
be independently controlled have been studied. These parameters are defined below . 

(i) The axial speed of the beam, expressed as a speed parameter c = c/{ca)\ where 
(ccr)i = v^Tq + TT^. It may be noted that c varies linearly with both Qp and M, as 
c = 2QpR/VTo + Tr^. 

(ii) The stiffness of the flexible pulley-support designated by the non-dimensional pa- 
rameter k. 

(iii) The amount of unbalance mei, generating the parametric excitation, expressed by 
the non-dimensional parameter eeo- 

(iv) The external damping parameter e/x. 

For the primary parametric excitation, Qp should be in the neighbourhood of 
which in turn depends on various system parameters, like c and k. The results are, 
therefore, obtained for the ranges of values of c and k, that yield — Qp- Further, 
results are presented for n = 1, i.e., when the first mode is primarily excited through 
the parametric excitation. 

Figure 5.2 shows the variation of the limit-cycle amplitudes with the axial speed 
for two different values of eeo- It is indicated, as discussed in Section 5.3.1, that when 
two limit cycles exist, the trivial one is unstable. On the other hand, when three 
limit cycles exist, the trivial one and the one with the highest amplitude are stable. 
The amplitudes of the non-trivial and stable limit cycles are seen to increase with 
increasing axial speed. Moreover, if three limit cycles exist, then a jump phenomenon 
with decreasing axial speed occurs, as exhibited in Figure 5.2. This feature has also 
been exhibited by a parametrically excited Duffing oscillator, as reported in reference 
[145]. Similar observations were also made by varying the stiffness parameter k. 

Figure 5.3 shows the effect of variation of the unbalance parameter eeo for two 
different axial speeds. At the lower speed, for which is negative, only the trivial 
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Figure 5.2: Variation of the limit-cycle £ 
o :eeo = 0.8, □ ceeo = 0.1; — stable, - - 



1 




limit cycle exists for low values of ceo- With increasing eeo, beyond a critical value, 
two limit cycles appear. When the axial speed is increased beyond a critical value, 
cr„ becomes positive. In such a situation, one, two or three limit cycles may exist 
depending on the value of eeo. At such high speeds, the possibility of a two-way jump 
phenomenon, as eeo changes, is clearly indicated in Figure 5.3. 

Figure 5.4 shows the effect of the external damping parameter at three different 
axial speeds. It may be noted, that with increasing damping, the non-trivial limit 
cycle vanishes at all speeds. The jump phenomenon is also exhibited as the damping 
parameter changes. Both one-way and two-way jumps are possible depending on the 
axial speed. 

The effects of the non-linearity on the stability of a parametrically excited axially 
moving beam are summarized in Figure 5.5. The parameters used in this figure are 
defined in Sections 5.2 and 5.3. Using equation (5.48) one can draw the curve AB 
delineating the stable and unstable zones as predicted by the linear analysis, where the 
zone marked II refers to the unstable zone. When the non-linearities are included, in the 
zone marked I, the equilibrium position turns out to be stable, whereas in zone II the 
equilibrium position becomes unstable and a stable limit cycle (with finite amplitude) 
is approached. However, in the zone marked III, the system approaches either the 
equilibrium position or a stable limit cycle depending on the initial disturbance. 

5.4 Near-Resonance Response with Simultaneous 
Parametric and Harmonic Excitations 

In this section, the response w{x, r) is obtained when an external harmonic force is 
present over and above the parametric excitation. Due to the simultaneous presence 
of both forms of excitations, the response usually turns out to be aperiodic. It is to be 
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Figure 5.3: Variation of the limit-cycle amplitudes with 'eeo- R 
0.05; o :cf = 0.89 , □ :c' = 0.94 ; stable, unstable. 








Figure 5.5: Stability boundaries of the parametrically excited system as predicted by 
the linear analysis. 
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noted that different combinations of frequencies of parametric and harmonic excitations 
are possible. These, in turn, may give rise to various complex patterns of the response. 
In what follows we shall consider a special case, where the exciting frequency is close to 
the natural frequency associated with the mode excited by the parametric excitation. 
Hence, along with equation (5.28) one can also assume that the forcing frequency Q,f 
is given by 

Qf=u}l + ea. (5.53) 

When an ^ (T, the response is almost periodic. In the absence of any parametric 
excitation, the steady-state response is also harmonic having the same frequency as that 
of the excitation. But in the presence of a small parametric excitation, the response 
can be assumed as 

W (x, r) = 2 -f- c.c. . (5.54) 

n=l ^ 

In view of equations (5.27) and (5.53), the following assumptions are justified : 


o>m — 0(e) = ea'jn (say), for m^n, 


(5.55) 


and 


with 


o-m = 0(e) = ea'^ (say), for m = 1 , 2 , 3 , ■ • • , 


(5.56) 


^an(r) = 0(e). ( 5 . 57 ) 

First, equations (5.54)-(5.57) are substituted into equation (5.25) and the harmonics 
are balanced (i.e., the coefficients ofe^r and e-W- are equated, separately, from both 

sides of the equation). Then using equations (5.30) and (5.31) and neglecting terms 
O(e^), one gets 






Lm^tn 


m=l 
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where 


= eNf + fi, (5_5g) 

Nr = {eo(l - ( 5 . 59 ) 

and 

fi = |Me'«/,o| . (5.60) 

Using equation (5.57) and orthogonality relations between the linear modes (see equa- 
tions (2.27) and (2.28)) in equation (5.58), and excluding terms O(e^), one gets 

^ + iOniQf - w„) = ^ $^Nf dx + ^ da:. (5.61) 

Further, the coefficients and o^jj’s can be obtained as, 


2J'oC(fi + rN,)(to 

- ^L)im 


if m^n, 


(5.62) 


and 


€S„ — a 


'm — -/rs 1 . .1 \ j . 


for m = 1, 2, 


(5.63) 


i(Q,f + 

Since Nf contains a term having explicit time dependence, equation (5.61) does not 
exhibit any periodic solution unless Cn — Let us consider the special case when 
= or, i.e., the frequencies of the parametric and external excitations are related. 
In this case, a steady-state solution to equation (5.61) exists. This solution can be 
obtained as explained below. 

Substitution of ^ = 0 in equation (5.61) yields the complex algebraic equation 

= + (5.64) 


For solving equation (5.64) to determine a„, one substitutes On = in the above 
equation and equate the real and imaginary parts from both sides to obtain 

Uni^f cos 6n = [G^ COS -f- sm On] 
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(5.65) 


—€■11 On sin ^ \Q^ + 
and 

an(^f -^n~ ^^n'S') sin 6n = ee^fin [C?^ cos 9n - sin 

+€(1 an cos dn + p [Q^ + co^M^] , (5.66) 

where 

S = 0t^/al, Q = Q^ + iQr = -col t da; 

Jo 

and 

M = I - Ko-i. / dr| + I E dx} . 

Now eliminating from equations (5,65) and (5.66) one gets, after some algebraic 
manipulations, a quintic polynomial in aj. This polynomial for a concentrated external 
force (at a; = xo), i.e., when f{x) is written as Fo6{x - xq), is given in Appendix C. 
For this particular loading, numerical results reveal one, three or five (positive) values 
of o2. However, a subsequent stability analysis only can confirm whether the values so 
obtained are physically realisable (i.e., stable) or not. 

The stability analysis can be carried out using the method outlined in Section 
5.3.1. However, an alternative procedure, considering an arbitrary perturbation to the 
solution, which gives rise to Mathieu type equation, has been carried out. 

5.4.1 Stability Analysis of the Steady-State Solution 

The stability of the steady-state solution, u;(a;,r) appearing in equation (5.54), is 
studied by perturbing it in the form 

Wp{x, r) = w(x, r) -f ws(x, r). ( 5 . 67 ) 
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Putting equation (5.67) into equation (5.25) and neglecting terms O(e^), one gets 


dWs 


where Wg = 


Nii = |2eo(l - k) cos 2%t^~ - 


(5.68) 


(5.69) 


( 6 .. 0 ) 

To analyse the stability of On, when -f the following perturbation is 

considered ; 


Wg - $„(x)^„(r) + C.C. . 

Substituting the value of w{x,t) from equations (5.54)-(5.56) one gets 

+ B$n^n + c.c = cLi + eL 2 + eLj-^ 4- c.c., 
dr dr 


(5.71) 


(5.72) 


where 


L, = {l,(s:,T),oy, t) = 1,2,3), 








ih. + 2 ( /' dx) 

da;2 Wo dx dx / dx^ J 




yo dx dx I dx2 I 




U. Kf/ji 

Z 2 (x, r) = 2eo(l — A:) cos 2QpT 
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and 


kix, r) = -lX(j)n. 

Applying the orthogonality relations (i.e., equations (2.27) and (2.28)) one finally gets 
the following equation and its complex conjugate: 


dr 


= i if! + i ilo 


dr 


d^n 


Assuming t' = fi,r + d„ and ri = S’ equation (5.73) can be written as 


(5.73) 


Qf 


dr' 


■ ^ = i + f (£ § 

t(/>»(^‘+C,)d.)?. + ^(/CE'sd.)§, , (5.74) 


where 


Lj = {lj(x,T'),0f y = l,2,3,.'-. 
Using equation (5.53), one can write 


T} = l + €rji-j . 

(5.75) 

Now, is assumed in the following series form: 



(5.76) 

Substituting equations (5.75) and (5.76) into equation (5.74) and equating the coeffi- 
cients of the like powers of e from both sides, one gets: 


(5.77) 


Jo ^ -^J dr' 
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The following solution of equation (5.77), 


(5.78) 


e>=6e‘n 

is first substituted into equation (5.78)., Then, to avoid secular terms in {»), the 
following condition must hold good : 

i% - r,)5 - TrS + Ge-^-b + i ‘ = 0 (5.79) 


where 


■ :« 1“ (r-S'i-*-) (/;&»■*) •• tf St -) tf ©>.*)i 


= ji i(r (/; a»-*) • (/;§ t*) n: /-■ 


For the existence of a non-trivial solution of 6 (= 6 * + iby) from equation (5.79) one 
finally obtains, (after some algebraic manipulations) 


771 = Fx ± [r^ - 2 r 2 + f-rl 

\</ \^nj ^ 


(5.80) 


with Ax = cos 29n + sin 20„. In the o„ - fi/ plane, equation (5.80) implies two 


curves 


n,« = wi(i-t>)h.j = i.2. 


where 77 ^^ corresponds to either of the values of Tji from equation (5.80) and denotes 

the associated stability boundary. 

For a fixed value of a„ and 0„, are calculated first < ilj], say) and then 
checked with the forcing frequency fi/. The amplitude is said to be unstable if 


fl/4 < 0 .f < fi/i* X 


(5.81) 
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and stable otherwise. 

It is to be noted from equation (5.80) that for some values of the parameters, rji 
does not exist. This implies a stable condition under such situations. 

5.4.2 Numerical Results and Discussion 

In this section, numerical results are presented to study the effects of several parameters 
on the steady-state response. All the results are obtained for To — 1, rco = 1/3 and 
Q/ ~ fip ~ aj[. 

At this stage it should be noted from equations (5.54)-(5.56) that the steady-state 
response of the beam is obtained (up to o(l)) as 

w(x,t) =ai<f>l{x) cos {QfT + 01 + pi), 

where tan pi = <j>l = The values of di and 9\ are now 

obtained by solving equations (5.65) and (5.66). Subsequently the effects of different 
parameters on the value of di are studied. 

Figures 5. 6-5. 8 show the response variable Si with increasing speed for two different 
sets of values of Fq and eeo. It is seen from Figure 5.6 that if the external forcing pre- 
dominates, then only one stable response is resulted. The value of Si, while remaining 
insensitive to the amount of unbalance, increases with increasing speed. On the other 
hand, if the parametric excitation predominates, then depending on the value of eeo and 
c , one or three unstable roots for Si may be obtained (Figures. 5.7 and 5.8). However, 
the usual jump phenomenon, observed in a harmonically forced Duffing oscillator, is 
not altered. Figure 5.9 clearly indicates that, if a flexible pulley-support is used, then 
the steady state response can be controlled by adjusting the phase relationship between 
the external and parametric excitations. 
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U'— ^ 1 L- "^1 I 

0.85 0.87 0.89 0.91 0.93 0.95 


c' 


Figure 5.7: Variation of the roots of equations (5.65) and (5.66) with d. R - 0.2, 
k = 0.5, 6f = 0.0, e/i = 0.05, Fo\/e = 0.1; o :eeo = 0.1 ; stable, unstable. 
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gure 5.8: Variation of the roots of equations (5.65) and (5.66) with c. 

= 0.5, $f = 0.0, €fi = 0.05, Fo^/e = 0.1; o -.eeo = 0.8 ; stable, - unstab 






Chapter 6 


STABILITY OF AN 
ACCELERATING BEAM 


6.1 Introduction 

The studies on .the non-linear vibrations of axially moving beams, carried out in Chap- 
ters 2 to 5, have been restricted to systems with uniform axial speed. Although the 
operating speed is normally maintained constant, the axial speed does vary consider- 
ably during the starting and stopping phases. It is the objective of the present chapter 
to bring out the effects of non-linearity on such an accelerating or decelerating beam. 

For a travelling string, numerical studies reveal that the linear response builds 
up or decays as the string is accelerated or decelerated, respectively [149, 150]. An 
experiment with a pipe carrying fluid was also reported in reference [150] to validate the 
numerical results. However, the stability of a beam travelling with constant acceleration 
or deceleration has received less attention, though the parametric excitation caused by 

periodic variation of the axial speed has been studied [70]. 

The equation of motion for non-uniform axial speed is given m Section 6.2. In 
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Section 6.3, the Lyapunov’s method is used to determine the sufficiency condition for 
asymptotic stability of a beam having any arbitrary acceleration. However, in some 
cases the method fails to ascertain the stability. An alternative method, based on 
Multiple Time Scale (MTS) is discussed in Section 6.4. Here, the change in the speed 
due to a constant acceleration or deceleration has been treated as a perturbation to 
the uniform speed.. Some numerical results are discussed in Section 6.5. 


6.2 Equation of Motion 


With the assumptions mentioned in Chapter 2, the non-linear equation of motion of a 
non-uniformly moving beam can be written in the following non-dimensional form: 


d‘^w d^w -2 m \ 
+ 2c^-^ + (c2-ro) 


d'^w d'^w dcdw _ 

dx"^ dx^ dr dx 




dx 


d'^w 
dx^ ' 


( 6 . 1 ) 


dr^ dxdr 

The non-dimensional symbols used in equation (6.1) have already been explained in 
Chapter 2. The boundary conditions due to simple supports are given by equations 
(2.10) and (2.11). 

When a small viscous damping term 5^ is present, the equation of motion turns 


out to be 

d^w -dw d^w 
9r ^ dxdr 




d'^w d^w dc dw 


dx'^ ' dx^ dr dx 
which can be recast in the standard state-space form as 

dW 

A-;^-|-BW-l-N = 0 
dr 

where A, B and W are given by equation (2.16) and 




d^w 


, -dw dcdw 


' (dW\‘ 

Jo [dxj 


dx 


d‘^W 

dx^ 


,0 


( 6 . 2 ) 


(6.3) 


(6.4) 


In the following, the stability of the response w{x, r) is analysed using either equation 
(6.2) or (6.3) whichever is convenient depending on the method used. 


136 



6.3 Stability Analysis Using Lyapunov’s Method 


In this section, Lyapunov’s direct method [151] is used to find the asymptotic nature 
of response of the system expressed by equation (6.2). 


The following is taken as a Lyapunov function , 



Assuming that the axial speed always remains positive (i.e., c> 0) and never crosses the 
first critical speed y/To + the chosen Lyapunov function is positive definite, since 
/o dx > Sq dx for a continuous differentiable function with iy(0) = 

u;(l) = 0, [152]. 

Equation (6.5) is differentiated along the curve satisfying the equation of mo- 
tion. Subsequent integration with the help of non-dimensionalized boundary conditions 
yields, 


dr 


\5r y 


i: 


dx — 2c 


yi dw dc / dw\ ^ ^ 

Jo dr dxdr ^ ^drJo \9xy dr Jo dx 


dc dw dw 
dr 


dx. 


It can be seen by simple integration by parts that 

ri dw d^w 


r 

Jo 


dr dxdr 

for the given boundary conditions. Thus, 


dx = 0, 


dV 

dr 


dc dw dw 
dr 


r^fdw\^ dc (dw\ , dc dw 

= -'/„ (^) ^-“d;/oUj 

Rearranging equation (6.6) as 

Jo I °'dT 2y/5 dr dx 


dx. 


dr 


f dc 1 [dcV] yi ( 

^^"^l^dr 4(5Vdrj j Jo \ 


dx , 


( 6 . 6 ) 


dx 


(6.7) 
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it can be seen that for a continuously accelerating beam, i.e., (dc/dr) > 0, 


— < 0 if 
dr 


Mc'\^ . /dc' 


(6.8) 


When (dc/dr) = 0, it is observed from equation (6.7) that (dl^/dr) < 0, i.e., the 
system is asymptotically stable. When (dc/ dr) ^ 0, the system is stable if 


(^) < ( 6 - 9 ) 

J_/ 

Further, it is seen that the beam is stable (i.e., ^ < 0) when both c < 0 and 
(dc/dr) < 0 and the inequality (6.8) is satisfied. The situation corresponds to an 
acceleration in the direction of -x. Thus, considering the symmetric nature of the 
boundary conditions, these two cases are physically identical. 

Although the stability is confirmed if condition (6.9) is satisfied, nothing can be said 
if the same condition is violated. Similarly, for a decelerating beam (i.e., (dc/dr) < 0), 
equation (6.7) can not be used to judge the asymptotic stability. This necessitates the 
choice of other functions. 


6.4 Stability Analysis Using MTS Method 


In this section, the stability of a uniformly accelerating or decelerating beam is analysed 
by the multiple-time-scale method. The change in speed is treated as a small pertur- 
bation to the uniform speed. Furthermore, the damping is also assumed to be small. 

Consequently, the speed of the beam and the damping factor are written, respectively, 
as 


c = Co -f ear (6.10) 

and 

and b - e^o. (6.11) 
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where |a| governs the magnitude of constant acceleration/deceleration. Using two time 

scales to = r and h = £t, and neglecting terms 0(e=), the partial differentials in time 
(r) can be written as 

±_d d 

dr ( 6 - 12 ) 

i 

(e.is) 

and the displacement as 


w{x,to,ti) - Woixytoyti) + ewi{x,to,ti). (6.14) 

Substituting equations (6.10)-(6.14) into equation (6.2) and equating the coefficients 
of the like powers of c from both sides, the following equations are obtained in the 
state-space form: 

t" ;A^+BM'„ = 0, (6.15) 

e' :A-^ + Bl^i-f-No = 0, (6.16) 

where . 


•NT _ /x ^^0 , d'^WQ , ^ d^Wo , ^ d'^wo , ^ ^ d’^wo , dwo 

“ 1 ''° 0(7 ^ + “■& 


fi¥)' 

Jo \ ox } 


^0 . 


Equation (6.15) can be solved as 


Wo(x,to,ti) = E 


2 ^ 


(6.17) 


(6.18) 


where anih) is the complex amplitude and the bar at the top denotes the complex 
conjugate. To solve equation (6.16), equation (6.18) is first substituted into equation 
(6.17) and then No is expressed as 


No = E Pn(x)e^"‘° + Q(a:,to) + c.c., 
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where c.c. denotes the complex conjugate of the previous terms and Q contains terms 
having all other frequencies except One readily finds Pn as 


f 1 CLn , . 

Pn(rr) = + + C0— ^ 

. , , dOn TT d^0n . . 

dir + “Td^ 


+ 


an“n 


('/« 


1 d.^n d <^^ ^ dVn 

0 drc dx / dx^ \Jq V 


dx 


dVn 

dx^ 


,0 


(6.19) 


Using the orthogonality relations between the linear modal vectors (see equations 
(2.27) and (2.28)) Pn is decomposed in terms of the linear modes as 


Pn(^) — [pmn^m(3i) ”1" 9mn^n(®)] , 


where 

/oCAPndx 

and 

/o‘CAP„<to 

;„‘?^A4„di' 

Now, the secular terms in the solution of equation (6.16) can be avoided if, 


Pnn — 0) 


i.e., 


r$S;^APndx = 0. (6.20) 

JQ 

Equation (6.20) can be expanded with the help of equation (6.19) as, 

- 2iCo [ ^5. d:.) ^ = 

A A J . TT /-Id^n-r , 

~ <Pn(pn<ix -TraOn — dx + ^CoQ!a„-^ / -3-5-0„dx 

^ Jo Jo dx wi, Jo dx^ 
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( 6 . 21 ) 


+ ia 




2 Jo dx 


4 d 




0 dx2 


dx . 


It is to be noted that for a travelling beam, 


lo =0, Im ^ ^„(^„dx =0 and Im =0. 

Now substituting a„ = o„e‘'’'‘ into equation (6.21) and equating, separately, the real 
and imaginary parts from both sides, one obtains 


i {2u,i£<l,J^dx-2icj^^<t.„dx'^^ = 


i I' - 2<Co I' ^ &) ^a„ = mocK jl ^?„dx 


( 6 . 22 ) 


. TT yl d^^n-T j •*■“3 /n 


dx dx ; 
r^-7 


S dVn 




(6.23) 


Writing /J dx = (5i > 0), equation (6.22) takes the following form: 

-wiycJo^^^n^ndx-ay^i. (6-24) 

Obvious from equation (6.24) is the fact that for a > 0, i.e., during acceleration, 
(da„/dtx) < 0 or the response decays. But for a decelerating beam (i.e., a < 0), the 
response may be stable or unstable according as 

Jo 


(6.25) 


141 



or 

f '$n<Pndx < (6.26) 

J 0 

respectively. However, the instability is only for a short time. The ephemerally ob- 
served characteristics of the instability can be explained as follows. As the beam 
becomes unstable, the amplitude grows but the speed also decreases since the beam 
decelerates. The numerically obtained value of 5i is found to decrease monotonically 
with decreasing speed. Consequently, with decreasing speed the beam eventually re- 
gains its stability as (da„/dti) becomes negative. Thus, the amplitude builds up to 
a limiting value after which it starts decreasing. This phenomenon was observed in a 
pipe carrying fluid when the flow was stopped [150]. It should be pointed out that the 
non-linear terms do not affect the stability but only change the frequency of oscillation. 
This can easily be seen from equation (6.23). 

6.5 Numerical Results and Discussion 

Numerical results are presented in this section for a beam having an initial tension 
To = 1. For a beam with its axial speed given by equation (6.10), the asymptotic 
stability is confirmed using the Lyapunov’s method (see equation (6.7)) if 

r a 
(Jo > 7-. 

4co 

The MTS method on the other hand predicts the stability for all possible values of 
So (So > 0). For a decelerating beam, the chosen Lyapunov function (see equation 
(6.5)) fails to ascertain the stability. However, it is evident from the MTS method 
(see equations (6.25) and (6.26)) that the stability depends on the values of So and 
|a|. Furthermore, it is obvious from equations (6.25) and (6.26) that in the parameter 
space (So, |a|) the boundary delineating the stable and unstable regions is a straight 
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line. Figure 6.1 shows these boundaries for three different values of co. The region 
above the line representing the boundary is stable. As expected, the danrping required 
to prevent instability increases with increasing initial velocity. 

The asymptotic or long-term stability, discussed in the previous section, is shown 
in Figure 6.2. The amplitude is approximately calculated as 

where 

A = ^ dx + aSi)/ dx - 2ico da:^ . 

When the damping is not sufficiently large, the response amplitude shows a temporal 
rise. The results have been obtained by assuming that only the first mode is excited. 
For the values of damping and axial speed considered, other modes are found to be 
stable. 

The effects of magnitude of deceleration on the maximum amplitude rise and the 
overshoot-time (i.e., the interval during whidi the response amplitude remains higher 
than the initial value) have been plotted in Figures 6.3 and 6.4, respectively. 
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|a| 


;y boundaries for a uniformly decelerating beam. 
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aximum overshoot of S', (in %) 



W\ 

Figure 6.3: Variation of the maximum overshoot of response amplitude with the mag- 
nitude of uniform deceleration. = 0.7, 
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Overshoot time 



Figure 6.4: Variation of the overshoot time with the magnitude of uniform deceleration. 
6o = 0.7. 
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Chapter 7 


VIBRATION OF A TRAVELLING 
BEAM HAVING AN 
INTERMEDIATE GUIDE 


7.1 Introduction 

The effects of the non-linear term on the free and forced transverse vibrations of a 
travelling beam under various operating conditions have been studied in the foregoing 
chapters. This chapter addresses to the problem of controlling these unwanted oscil- 
lations. A simple passive controller, in the form of an intermediate guide has been 
suggested. To analyse the effectiveness of such a controller, the response of a trav- 
elling beam to a harmonic excitation is considered with and without the guide. The 
equation of motion is derived in Section 7.2. The complex normal modes of the con- 
strained system are obtained in Section 7.3 by neglecting the non-linear term. The 
near-resonance response is then obtained through a modal analysis. Both the stiffness 
of and dry friction in the guide are taken into consideration. Since the frictional force 
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frictionlcss guide 



Figure 7.1: Schematic diagram of a travelling beam with an intermediate guide. 

between the guide and the beam has a destabilizing effect, rollers can be used to main- 
tain the contact between them. Similar analysis is carried out in Section 7.4 including 
the non-linear term. Numerical results are included to highlight the salient features. 

7.2 Equation of Motion 

Consider a slender beam moving axially between two frictionless guides as shown in 
Figure 7.1. The axial speed, c*, is maintained by means of two rigidly mounted pul- 
leys. An intermediate guide, consisting of two free rollers, with finite compliance is 
held in contact with the beam by means of an external force JV*/2, provided by the 
precompression of the compliant member. In addition, as indicated in the figure, the 
stiffness of the guides is assumed to be same in both the directions. 
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Under the usual assumption of small longitudinal vibration (in comparison to the 
transverse vibration, i.e., assuming u* = 0{w*^)), the non-linear equations of motion 
for the coupled vibration in the axial and transverse directions are 

dw* d^w* 

(7.1) 
and 


pA 


d^u 


dt^ 


* * d^u 

+ 2c* 


+ d 


id^u* 


d^dt ■ de \ 


pA 


d^w* 


at2 


+ 2d 


d^dt 


+ c* 


td^w* 


rputi 

° de 


d^w* _ _ d^w* 
+ EI,- 




a^2 j 


d .du* dw* 

(7.2) 


where JF\* and Fj* denote the frictional forces and are given by 


FI 


= Pf 


F* = pj 





(7.3) 

(7.4) 


In equations (7.1)-(7.4) , Fjf/2 and pf denote the stiffness of the guide placed at ^ = d*, 
and the coefiBcient of dynamic friction, respectively. The other symbols are explained 
in Chapter 2. 

As discussed in Chapter 5, the vibration coupling between the two spans is ne- 
glected. The boundary conditions can then be written sis 


u*{0,t) = u*{l,t) = Q, 

(7.5) 

w*{0,t) = w*{l,t) = 0 

(7.6) 

aV( 0 ,t) d!^w*{l,t) 

9^2 - 0- 

(7.7) 


Using the following non-dimensional parameters, 

Fi = F-KEA-fY, i = 1, 2, N = N'KEAfY K, = KyUEA), d^d-/l, 
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along with those given in Chapter 2 (see equation (2.6)), equations ( 7 . 1 )-( 7 . 7 ) can be 
written, respectively, as 


5r2 ^ dxdr dx^ 


1 0/1 om .dw d^w 


d'^w 

^ dxdr 


+ (c^-To)— + 


dx^ ' dx^ 

—KfwS{x — d), 

^ \N Kf , 
=M/ [j + -^‘^id,T) 


I dudw .dw.^ 
7 ^ dx dx 2 9x 


F 2 = jif 


N Kf 


w{d, t) 


L2 2 

ti( 0 ,r) = u(l,r) = 0 , 
u;(0,t) = w{I,t) = 0 


(7.9) 

(7.10) 

(7.11) 

(7.12) 

(7.13) 


and 


(7.14) 


d^w(0,T) _ d^w{l,T) _ 
dx^ dx^ 

For a small value of 7 (i.e., 7 < I), the longitudinal inertia in equation (7.8) is neglected 
to yield 


I d^n ndwd^w 


- {Fi + F2)6{x - d). 


72 5 a ;2 ' dx dx'^ 

which, when integrated twice with respect to x, results in 

2 


u(x, r) = - Y (^) ^ ~ 


where 


H(x-d) = 0 ; x<d, 
= I; x>d. 


Using equation (7. 1 2) 


the unknown constants of integration /((r) and /^(r) are ob- 


tained as 




(7.16) 
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and 


= (7,17) 

Now combining equations (7.10), (7.11), (7.15)-(7.17) and (7.9), the equation of motion 
for the transverse vibration can be written in the following simplified form: 


+ 2c^ + - {T. - - <0 + - ci)}] ?? + 




dx^ dx‘^ 


,dw 


~ — d) — e 




d'^w 


(7.18) 


2 

where €(= ^) is a small parameter, i.e., c < 1. This equation can also be written in 
the familiar state-space form as 


Ai 


dW 

dr 


-t-BiW^ = €N, 


(7.19) 


where Ai 


1 

0 


G Ki ' 


dw 


,Bi = 


,W=< 

dr 1 

1 

OKi 


1 

1 

0 

1 


1 ^ . 


V, Ki = (c^--Tb)^ + | 


a* j a* 


fijNH{x d)-^ — 6{x — d){fj,fN^ - Kf), G = 2c^, I is the identity operator and 


N 


' 

Ja \ 5 a: 


da: 


d^w 

dx^ 


,0 


(7.20) 


It is not difficult to verify that Ki is a self-adjoint and positive definite operator. 
This fact is essential to obtain the orthogonality relations between the complex normal 
modes as discussed in the next section. 

In order to facilitate computation, the non-analytic functions involved in equation 
(7.18) can be eliminated by writing the equation separately for two domains, namely, 
0 < a:i < d and d < 0:2 < 1 as 




+ 2c' 


dxjdr 


+ [c" - Tj] 


d'^Wj 

dxj 


d^Wj 

dx] 



(7.21) 


152 



for j = 1,2, where Wj is the transverse displacement in the j-th. domain and 


Ti = To-^fN{l-d), 

T2 = To + HfNd 


or 

AT = T2-Ti = HfN. 


The matching conditions aA x = d can now be written as 


Wl(d,T)=W,(d,T), 

and ^(d,T)-^{d,T) + Ar^(d,T)-ir/<»iWr) = 0. 


(7.22) 


The non-linear response of the beam can be studied by solving equations (7.21) and 
(7.22). Since the non-linear term can be taken as a small perturbation to the linear 
equation of motion, the response of a linear beam is discussed first before taking up 
the non-linear one. 


7.3 Free and Forced Responses of the Linear Sys- 
tem 

In this section, the free and forced responses of a linear travelling beam (i.e., with e 0 
in equation (7.19)) are discussed. As already mentioned in Chapter 2, no stationary 
mode shape exists for a travelling beam, but the harmonic oscillation is still possib 
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some fr6qu6nciGS, known ns th6 ^nnturnl frnqucncics . Considering the complex normnl 
modes, the natural frequencies and the complex mode shapes can be obtained. The 
response of the beam, w{x, r), at any one of the natural frequencies can be written 

as 

. /7 . f 

(7.23) 




W{x,t) = f 


with $(x) = 


^ J 


, where (f>{x) is the complex normal mode shape. Substituting 


equation (7.23) into equation (7.19) (with e = 0) and equating the coefficients of e" ' 
and separately from both sides, one gets 


iui^Ai^ + Bi$ = 0 


(7.24) 


and its complex conjugate, respectively. It is to be pointed out that the partial deriva- 
tives appearing in Ai and Bi (see equation (7.19)) are. replaced by total derivatives in 
equation (7.24) which is now solved together with the boundary conditions to deter- 
mine w^’s and (j>’s. The n-th natural frequency and mode shape so obtained are denoted 
by o;^ and (f)n, respectively. 

It should be noted that the above solution is to be obtained numerically. In stead 
of solving equation (7.24), which contains several non-analytic functions, it is broken 
in two domains as explained in Section 7.2. Thus, equation (7.24) is reduced to 


-(J)% + 2icw‘^ + (c» - + 


dV; 


dxj- 




dxj 


0 , 


(7.25) 


with j—1, 2. The matching conditions at x = d are still given by equation (7,22) with 
the partial derivatives replaced by the total derivatives and w replaced by 4>. Assuming 
the solution in the form 


(«i) = E ajfceP-#**'' , i = 1, 2 

fc=i 
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and applying the boundary conditions, one obtains 


( ^ 

Oin 

> = [ATil < 

ai3 

> 

. “12 , 


014 ^ 

[ Q!21 

1 

t = m . 

^23 

► 

( a22 J 


h ^24 , 


[i^o] 


= 0 . 


(7.26) 


and 


The four unknowns ais, an, 0:23 and 0:24 are now substituted in the matching conditions 
to obtain a relation like 

f 

«13 
Q!14 

Q!23 
0!24 

The existence of a non-trivial solution implies 

det [JiTo] = 0- 

Since the matrix [/To] is a complex one, both the real and imaginary parts of the 
determinant must vanish simultaneously. The values of w' and pjkS (j, k = 1,2, 3, 4) 
are obtaind numerically by solving equations (7.25) and (7.27). 

Considering the symmetric and antisymmetric nature of the matrices Ai and Bi, 
respectively, the following orthogonality relations amongst the complex modes and 
their conjugates are easily obtained: 

(7.28) 


(7.27) 


r da: = 0 for all m and 

Jo 


n. 


and 


j' = 0 

The above normal modes are used, as m Chapter 2, to derive the response otsuA 
a beam, when an external force f{x,T] is applied. Following the procedure detarled 
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in Chapter 2, the steady-state response of the beam to a harmonic excitation, i.e., 
f{x, t) = f{x) cos CIt, can be obtained by solving the following equation of motion: 



Ai^ + BiW = f, 

(7.30) 


with f = {/(a:,r), 0}^. 

(7.31) 

The steady-state response W (x, r) is finally given by 



Wix,r) = ^ f;(Pn^n + C.C., 

^ n=l 

(7.32) 

where 

. ig 2S}'^lfidx 

(7.33) 

and 


(7.34) 


with fi = {f{x)/2, 0}^. It is readily seen that at or near resonance the magnitude of 
Pn will be very high. Hence the aim of controlling the vibration is to reduce the value 
of Pn. 


7.3.1 Numerical Results and Discussion 

Numerical results obtained from the linear analysis (presented in the previous section) 
are now discussed. All the results, unless otherwise mentioned, apply to a travelling 
beam having a tension T 2 = 1. The important physical parameters which control the 
performance of the system are 

1) the speed of the beam, expressed as a parameter d = c/(c«:)i where (ccr)i = 

+ T2, 

2) the location of the guide, d, 
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3) the stiffness of the guide, Kf and 

4) the coefl&cient of friction between the beam and the guide, ixf. 

As shown in Chapter 2, an increase in the axial speed reduces the natural frequencies 
of vibration of a beam without any intermediate guide. The same effect is also observed 
even in the presence of a guide. In Figure 7.2, the first and second natural frequencies 
are plotted against the axial speed for various values of Kf. For the chosen guide 
location, i.e., d = 0.5, the second natural frequency, Wj, as expected, is insensitive to any 
change in the stiffness Kf. However, the first natural frequency and consequently 
the first critical speed considerably increases with increasing Kf. Thus, for large values 
of Kf (i.e., when the guide behaves like a rigid support), the first natural frequency of 
the guided beam tends towards the second natural frequency of the unguided beam. 

Figures 7.3 and 7.4 show the variation of the natural frequencies with the guide 
location. It is observed that the first natural frequency u[ is maximum with d = 0.5, 
whereas the value of Wj attains a maximum when d = 0.25 and d - 0.75. For a 
stationary simply-supported beam, the point a; = 0.5 corresponds to, the antinode of 
the first mode and the node of the second mode. Although the concepts of nodes 
. and antinodes can not be used for a travelling beam, the response envelope, however, 
attains a maximum or minimum at a: = 0.5 when the beam vibrates at its first or 
second natural frequency, respectively. Similarly, for the second natural frequency, the 
response envelope reaches a maximum value at a; = 0.25 and x = 0.75. Thus, from 
this example, it is confirmed that the guide makes the system most stiff if it is placed 
where the response envelope reaches its maximum. Attention may be drawn to one of 
the general results for a constrained flexible system, known as the ‘eigenvalue-inclusion 
principle’. For a gyroscopic system, it has been shown [30] that for a beam or string 
having a stiffness constraint, the n-th natural frequency of a constrained system (say 
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c' 


Figure 7.2: Variation of uj[ and w* with c . d = 0.5, AT = 0, : a;| with no guide; 

: iui, Kf = 10, : u}[, Kj = 100, — : u)\, Kf = 1000, — x— : ui, 

Kf = 0 , 10 , 100 , 1000 . 
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n^) satisfies the following inequalities 

where is the n-th natural frequency of the unconstrained system. Numerical 
results obtained for the guided travelling beam are in conformity with the above result. 

Figure 7.5 shows the variation of (jj[ with the coefficient of guide friction, fj,/. It 
is seen that the friction reduces the stiffness of the system. Prom equation (7.21), 
it can be concluded that the friction introduces a compressive load in the first span 
(0 < rr < d). This compressive load may become so high that the span may undergo 
divergence instability, as shown in Figure 7.5. Thus, in addition to damaging the beam 
surface, the presence of friction may cause instability in the system. 

For a guided travelling beam no ‘frequency-loci veering’ has been observed. This 
supportes the observation reported in reference [47]. The non-cyclic disordered beam 
discussed in this reference corresponds to the present model with Kf —¥ oo and d — 0.5. 
However, the phenomena of ‘frequency-loci veering’ and ‘mode localization have been 

observed for a guided travelling string [29]. 

. The effect of the intermediate guide on the steady-state harmonic response is shown 
in Figure 7.6. Only the magnitude of pi is plotted for a beam excited by a point 
harmonic load f{x) = FoS{x-xo) at xq = 1/3 with the excitation frequency Q « By 
changing the natural frequency, the guide helps to avoid the resonance at a particular 
speed. As expected, the near-resonance response decreases with increasing stiffness of 

the guide. 

Figures 7.7 and 7.8 show the effects of the guide-location on the steady-state re- 
sponse. As seen from these figures, the guide when suitably placed can attenuate the 
response level of a resonantly excited beam. As the application of the guide shifts the 
natural frequency, the resonance condition can be avoided. If the guide is flexible, then 
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its optimal location co^esponds to the position where the response envelope attains 
a maximum and depends on the mode which is resonantly excited. In- other words, 
to minimize the vibration, the guide is placed at a location where the response of the 
uncontrolled system is maximum. However, as the stiffness of the guide increases, the 
location of the guide looses its importance as long as it does not coincide with any of 
the nodal points of the response envelope. The shift of the natural frequency, after 
placing the stiff guide anywhere within two successive nodal points, is so large that the 
response amplitude becomes practically insensitive to the guide-location. 


7.4 Effects of Non-Linearity 


In this section, the effects of non-linearity on the free and near-resonance forced re- 
sponses are presented. It may be mentioned that the relationship between different 
linear natural frequencies depends on the various system parameters. Some typical 
relationships may give rise to internal and combinatorial resonances, which are not 
considered in this work. Thus, only such combinations of parameters are taken where 
all the natural frequencies are distinct and no special relation between them exists. For 
such a system, the free and near-resonance forced harmonic responses are analysed. 

In the absense of any internal resonance, the harmonic response of a free travelling 
beam has already been presented in Chapter 4. Following similar method first the 
complex non-linear normal modes are obtained for this constrained system. For free 
vibration in the n-th non-linear normal mode, the response is assumed as 

Wix, r) = + pr.ix)e-^’^\ (7.35) 


where 




>, with ipn and a;„ as the n-th non-linear mode shape and the 


corresponding frequency, respectively. Both these quantities, however, are amplitude 
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dependent. Following the method outlined in Chapter 4, a;„ 
respectively, as 

and can be derived, 


Wn = wi + e;d{”^ + ..-, 

(7.36) 


+ cAi H , 

(7.37) 

where 

pin) _ 

Ml UjUa , 

^ ‘'71 

(7.38) 


00 

Al = aa , 

mj^n 771=1 

(7.39) 


a z=z ^ ^n) , 

ym A 7 j~- ^ ^ 


h — ^ _ 1 ft Q 

^m A l \ i JL 5 ^ *“ ^9 "j ^5 

4a;^+a;^ <„ 

*n=/ ¥^Ai$„dx, n = l,2,3,--- 

JR 


and 

A(^m An) = 2 (j 4i 


‘ . dV~ 


0 d X 


d^n d^„ 

dx dx 


dx + 


dx . 


The above-mentioned non-linear complex normal modes can be used to obtain the 
near-resonance harmonic response of the beam when excited by an external force 
/(x) cos fir, i.e., to solve the equation of motion 


Ai^ + B,-tN = f. 

or 

where f = {/(x) cos Or , 0}^. 

The principal harmonic response of the beam can be assumed as 


(7^40) 


Wix, t)^ 1 {A{x)e^^ + A(x)e-‘''") 


(7.41) 
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For the near-resonance excitation « ujL the following response is obtained 


A = an'^n + ^ ^ ^ ^ 

m^n m=l 


(7.42) 


One can follow the steps detailed in Chapter 4 (see Section 4.3) to get the following 
results, valid up to 0 (e): 


OjTI ““ 




i(fl - a>„) * jAit„ dl’ 


(7.43) 


— 


2/o 


and 


ebm = 


»(n - ij„) si <k 

2;o‘47fjdc 


; m^n 


; m = 1, 2 


where fi = {f{x)/2, 0}^. The complex cubic equation (7.43) can be solved numerically. 
Assuming 


Qifi — 0,fiB 


equation (7.43), for a point load f(x) = F(i5{x — a:o)> takes the form of equation (4.35) 
which is reproduced below: 




- ^i) 


- i/ ((Mf )2 -h (M/)2) 


A^+ 


K)^(l-ri)^ 

52 


2u(M^Q^ 4- MiQ[) 


[(Qf)“ + («{)"] = 0 , 


(7.44) 


where 


= <^n/da: = Qf + i(5(, 

•^1 ~ 9m^m ^ I 53 ^^1 ~ 2 M/, 

/3{"^ _ f) _fFo\^e , _ 
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It may be mentioned that depending upon the system parameters, either one, or three 
roots may exist. As shown in Chapter 4, the intermediate root is always unstable 
and is not observable in reality. It was also shown in Chapter 4 that the steady-state 
response of the beam, up to the order o(l), is 

w{x, t) = an(f>l{x) cos {Q.T + 6n + Pn), 
where tan Pn = 4>i = \/ (^f 


7.4.1 Numerical Results and Discussion 


In this section, numerical results showing the effect of non-linearity on the free and 
forced responses of a guided travelling beam are presented. The tension of the beam, 
Ta, is again taken as unity, i.e., Ta = 1. 

Owing to the presence of the non-linear term, the natural frequency shows a hard- 
ening characteristics, since the term (see equation (7.36)) is positive. In Figure 

7.9, the variation of (aa with the guide location, d, is shown. As seen from the 
figure, the non-linear effects depend strongly on the location and the stiffness of the 


guide. 

As is well known, the Unear theory, when compared with the non-linear theory 
for such a hard system, over-estimates the steady-state response when the excitation 
frequency is equal to the natural frequency i.e,. fl = wi and under estimates the 
response when O. > u,‘„. Therefore, special attention should be given to the response 
for n > Figure 7.10 shows the variation of amplitude with the guide location. 
Where more than one value of the amplitude A are obtained from equation (7.44) 


(with n = 1), the higher stable amplitude is plotted. Three excitation frequencies 
n = 8.7, 12 and 15 are considered. These frequendes are such that the first mode is 
primarily exdted. The significance of the choice of the guide location, so far as the 
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maximum amplitude is concerned, can be clearly seen in Figure 7.10. It is seen that for 
Q = 8.7 (i.e., = 10.87), the non-linear response is less compared to the linear 

response. If a guide is placed, the linear natural frequency increases (see Section 7.3) 
and the excitation frequency becomes much less compared to the natural frequency, 
i.e., the frequency at which the resonance occurs. For a suitable guide-location, the 
difference between the excitation and the natural frequencies becomes so large that the 
effect of non-linearity is hardly perceptable. When (for example, f2 = 12 or 

15), the response amplitude predicted by the linear theory is much smaller than that 
obtained by the non-linear analysis. If a guide is now placed, u[ approaches and 
the response amplitude decreases which is a typical characteristic of a hard system. It 
is to be pointed out that even if the linear theory may not suggest the requirement 
of a guide by under-estimating the response, a guide may be required because of the 
presence of the non-linear term. 
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Chapter 8 


CONCLUSIONS 


8.1 Conclusions 

Different aspects of non-linear vibration of a travelling beam are studied in this thesis. 
The effects of the geometrical non-linear term taken as a small perturbation to the 
linear equation of motion of the beam are obtained by extending the results of lin- 
ear vibration. Both complex modal analysis and wave-propagation theory are used to 
study the free and forced vibrations of a linear travelling beam. An idea of ‘non-linear 
complex normal modes’, similar to the ‘non-linear normal modes’ for axially station- 
ary systems, has been introduced. These non-linear complex normal modes are used 
to obtain the near-resonance response under an external harmonic and/or parametric 
excitations. For a non-resonantly excited beam, the non-linear response is obtained 
using the results of the wave-propagation analysis. The effects of non-linearity on a 
non-uniformly moving beam has also been studied. A simple vibration controller in 
the form of a roller guide has been suggested. The major conclusions of the analyses 
presented in the thesis can be summarized as listed below: 

(i) The natures of the waves associated with the free modal vibration of a travelling 
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beam change depending upon the relations between various system parameters. During 
linear modal vibration, either two propagating and two evanescent waves (one of each 
kind in both upstream and downstream directions) or four propagating waves (one in 
the upstream direction and the remaining three in the downstream direction) appear 
depending on the axial speed. The well-known ‘phase-closure’ principle is satisfied for 
both the cases. The ‘phase-closure’ principle is also satisfied during non-linear modal 
oscillations of such a beam. 

(ii) The near-resonance response to an external harmonic and/or parametric excitation 
can be obtained easily using the non-linear complex normal modes. This method is 
computationally more eflacient than the commonly-used Galerkin’s technique. 

(iii) The steady-state response of the beam to a non-resonant hard excitation can be 
obtained by wave-propagation analysis. In this method, a closed form transfer function, 
i.e., the Laplace transform of a point-impulse response, is first obtained. The linear 
response obtained using this transfer function is then perturbed to get the effects of 
the small non-linear term. The present method is shown to be more accurate and 
computationally more efficient than either Galerkin’s technique or non-linear complex 
modal analysis. 

(iv) The non-linear natural frequencies of a travelling beam increase with the amplitude 
of vibration. The ‘jump’ of the amplitude of the response envelope can be observed un- 
der a near-resonant excitation. Further, one or multiple ‘limit cycle’s of the amplitude 
of the response envelope are possible when the beam is parametrically excited. These 
three phenomena show striking resemblance with those exhibited by a hard Duffing 
oscillator. 

(v) For a parametrically excited travelling beam with an external harmonic excitation, 
the effects of various system parameters depend on the relative strengths of the two 
forms of excitation. Moreover, a particular phase relationship between the two forms 
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of excitations (whose frequencies are related) results in a nainimum value of the steady- 
state response. 

(vi) A continuously accelerating beam remains stable, but for a decelerating beam in- 
stability may appear depending upon the magnitudes of damping and deceleration. 
However, the amplitude of a decelerating beam does not grow unboundedly and the 
beam regains its stability in the long run. The non-linear term does not play any role 
so far as the stability is concerned. It merely changes the frequencies of oscillations. 

(vii) An intermediate guide can be effectively used as a passive controller of vibration. 
The choice of the suitable guide-location plays a very important role in reducing the 
level of vibration. The guide also has a stabilizing effect so far as the divergence in- 
stability is concerned. But the frictional force between the guide and beam may add 
to instability. Thus, to use the controller effectively, the friction has to be minimised. 
Towards this end, a roller-guide is suggested. 

8.2 Scope of Future Work 

The foilwing are the directions along which future researches can be motivated; 

(i) The vibrations of the travelling beam having a single span have been considered 
in this thesis. A natural extension of these studies to multi-span beams (for example, 
band/wheel systems) can be carried out. 

(ii) The equilibrium configuratiou of the beam has been taken to be trivial. However, 
the case where the beam has initial transver^ dehection deserves further attenition. 

(iii) The study can be extended to the beam's having imperfections like weld melts etc.. 

(iv) Chaotic behaviour and the associated statistical properties of the non-linear re- 
sponse have not been considered in the present work. These can be take 
subjects for further investigations. 
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Appendix A 


Proof of Identity (2.23) and (2.47) 


The identity can be easily proved by first noting the following fact [153]. For a block 
matrix 


[A/l = 


wi [M2I 

[JWsl [M4] 


J 


det[M] = det[Mi]det[M 4 ]det (l - [Mi] ^[M 3 ][Mi]"^[M 2 ]) . 


Assuming the members of matrix [M] as, 




gtfc4 

, W! = 

^ikz 

1M3] = 

1 1 ' 
kl kl _ 

, [Ma] = 

1 1 

kl 4 _ 

) 


and performing a few row operations, one can easily verify that equation (2.23) turns 
out to be identical with the following : 


detfilf] = 0. 

Further, it can also be noticed that [i? 4 ][i 23 ][-R 2 ][J?i] = [Mi]-^[Mz][Mi]-^[M 2 ], whexe 
[RjYs are the same as appearing in equation (2.47). Thus, det[M] = 0 also implies 
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equation (2.23). 
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Appendix B 


Forced Response of a Travelling 
String 


For a linear travelling string, equation (2.60) is replaced by 

s^w + 2sc^ + (c^ - 1)^^ = 5{x- xo)e~‘‘^°. . (B.l) 

The response, from the wave-propagation analysis, formulated in Section 2.4.2, is ob- 
tained as 


ih(x,s) = fora:>a;o, (B.2) 

= -I- C 2 e^*® - fora;<xo, (B.3) 

where Ai = iki = 5/(1 - c), A 2 = ik 2 = -s/(l -I- c), Di = -e-^^°/2s, D 2 = 6-“^ /2s. 

Applying the boundary conditions t£)(0, s) = w{l, s) = 0, the following expressions 
of Cl and C 2 are obtained: 


Cl = _ eA2(i-xo)j ^2s{e^^ - 


and 


C2 = Je'^2(i-®o)+Ai _ gAi(i-xo)l f2s{e^^ - e^*), 
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which, when substituted into equations (B.2) and (B.3), yield 

w{x, s) = |^e'^2(s-a:o) _ gA2i— Ano _ gA2(l-a:o)-Ai(l-x) gA2-Ai(l+a:o-x)j 


-4T0 


X- 


o-sra 


X 


for a: > xq. 

(B.4) 

-Ai(xo-x) _ gA2(l-xo)-Ai(l-i)j 


for X < xo- 

(B.5) 


2s(l — e-^2-Ai) 

For To = 0, the above equations are identical with the transfer function obtained by 
Yang and Tan [60]. As mentioned therein, s = 0 is not a singular point since w{x, s) 
is finite as s 0. The only singular points are s = ±iUn = ±m7r(l — cP), when 
Ai = ±m7r(l + c) = and A 2 = — c) = A^”^ 

The temporal response is obtained by inverse Laplace transform of equations (B.4) 
and (B.5) as: 


w 


or 


{x,r) = Y. 


n=l 




\(«) 






2(1 - + 2s (^ - 5^) 


e‘““(T ^^+c.c. for T > To, 


r.A<’*)x _ gA[">xWg-A(">xo _ 1/1 /T, 

— y' ^ ig»W„{T-To) _|. g g £qj. ^ 0) 

^ ^ iimr 


00 /^Ali-'^x 

E 

n=l 

Since for a travelling string, the linear normal mode <l>n can be written as <f>n = 
gA["’x, Qj^g gg^jj verify that 




(B.7) 


H [ KK - 2fc dx = inir. 

Using equations (B.6) and (B.7) the response is finally obtained as 

/ \ 4>n{^)<l>ni^o) gtwj,(T-To) ^ g g_ fgj. (B.8) 

2u}\, fo^ Mn - 2ic /o^ ^ 

Equation (B.8) is identical to equation (2.55). 
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Appendix C 


Algebraic Equation for Obtaining 
the Near-Resonance Response 


The final algebraic equation, after eliminating from equations (5.66) and (5.67) for 
f{x) = Fo6{x — xq), is obtained as 






- u {(M^f + 


+ 


n 

52 


A A^ 


4(% - 
52 




j - 2v(Q^M^ + Q{M{) + ~ + {Ml)%} 


+ A‘ 




2v 


+ -^{hQiM^ A hQ{Ml) A - /5I3) 


+ 


S 

2v 


ik - k)iM,^Qi A Q^Ml)] -f- A 1 + |: {kiQfy 


+ hiQi?} - ^QfMf (12 + ll) _ |q{m/(Z2 + il) 

+ Y^iQiih - h) + + QfMDihh - hh) 


A 


+ ^5) + (Q()"(?^ + ll) - 2 Q^Q{{l 2 U - ^3^5)} 


= 0 , 
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(C.l) 


where Mi = M/Fq, Qi = Q/Fq, A = ml, h = [(G«)2 + (G^^J+eV'*, 

I 2 = (Qf - u}l) + ee'oG^, h = (% - - te^G^, k = ep' - ee'oG^ k = eix + ee^G^ , 

V - e{FQ/tnY/S'^. 
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